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1.0  INTRODUCTION 


Most  studies  of  cloud  motions  and  morphology  have  been  done  in  the 
Eulerian  frame,  that  is,  a  frame  fixed  with  respect  to  the  earth.  For  most  such 
studies,  the  Eulerian  frame  is  indeed  the  most  meaningful.  For  tropical 
cyclones,  however,  it  would  also  be  useful  to  study  cloud  motions  and  morphology 
in  the  Lagrangian  frame,  that  is,  a  frame  centered  on  the  storm  (and  therefore 
moving  with  it)  and  oriented  with  a  principal  axis  in  the  direction  of  motion  of 
the  storm.  Such  studies  would  hopefully  reveal  new  insights  into  the  behaviour  of 
tropical  cyclones. 

With  the  advent  of  cloud  imagery  from  the  Geostationary  Operational 
Environmental  Satellites  (GOES)  and  similiar  satellites,  this  type  of  study  is  now 
possible,  for  it  requires  a  long  series  of  images  taken  at  regular  intervals, 
preferably  intervals  as  short  as  possible.  With  GOES,  the  images  can  be 
generated  every  half  hour  for  images  of  the  entire  disc  of  the  earth,  or  every  15 
minutes  if  only  a  portion  of  the  full  disc  is  viewed.  From  these  images  subsets 
can  be  derived  that  are  centered  on  the  storm. 

The  purpose  of  the  contract  discussed  in  this  report  was  to  develop 
software  that  would  transform  tropical  cyclone  images  (from  GOES)  into 
Lagrangian  coordinates,  create  interpolated  images  to  replace  missing  or  garbled 
ones,  and  transfer  the  resulting  images  to  videotape  for  later  playback  through  a 
normal  TV  monitor.  The  images  provided  by  the  government  were  of  Hurricane 
Ignacio,  which  occurred  off  the  west  coast  of  Mexico  during  October  1979.  The 
images  provided  were  spaced  two  hours  apart  and  covered  a  6.5  day  period.  In 
all,  76  good  images  were  provided;  three  new  images  had  to  be  generated  by 
interpolation  to  replace  two  missing  and  one  garbled  image.  Thus,  a  total  of  79 
images  in  the  Lagrangian  frame  were  produced  and  transferred  to  videotape. 

In  the  next  section,  we  discuss  the  software  (including  procedures  for  using 
it)  and  production  of  the  videotape.  The  third  section  discusses  the  results. 
Finally,  Appendices  A,  B  and  C  present  flow  charts  and  listings  of  the  software. 
A  list  of  the  abbreviations  used  in  this  report  is  contained  in  the  final  appendix 
(D). 
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2.0  IMAGE  PROCESSING 


The  transformation  of  an  image  of  a  tropical  cyclone  into  Lagrangian 
coordinates  and  final  storage  onto  videotape  is  a  process  that  has  many  steps. 
The  process  involves  the  use  of  the  Data  General  Computers  associated  with 
SPADS  (either  the  Eclipse  S/250  or  the  Nova),  the  CDC  computers  of  FNOC 
(SPC  and,  possibly,  HAL),  and  a  special  piece  of  hardware  known  as  an  RBG  to 
NTSC  color  encoder.  The  entire  process  is  summarized  in  the  flow  chart  shown 
in  Figure  2-1.  The  following  sections  discuss  the  software  referred  to  in  this 
figure  as  well  as  the  procedures  for  using  it  and  for  transferring  the  images  to 
videotape. 

2.1  Initial  Processing  on  SPADS 

The  initial  processing  of  the  image  uses  software  previously  developed  for 
the  SPADS.  This  software  is  documented  elsewhere  and  will  only  be  briefly 
mentioned  here.  The  exception  to  this  is  the  program  ANGLEL,  which  was 
developed  specifically  for  this  contract;  it  is  discussed  toward  the  end  of  this 
section. 

The  user  must  start  with  a  large  sector  image  file  in  a  format  compatible 
with  the  SPADS.  This  image  file  can  either  be  directly  captured  from  GOES 
using  GADHS  or  it  can  be  taken  from  magnetic  tape  in  the  MCIDAS  format  and 
converted  to  the  format  acceptable  to  SPADS  using  the  program  MCODIS  or 
MCODIS2.  This  image  sector  is  then  subset  into  a  512  x  512  image  using 
GIMAGE;  the  starting  location  (upper  left  corner)  of  this  subset  is  chosen  such 
that  the  center  of  the  tropical  cyclone  is  close  to  the  center  of  the  image.  (It 
may  be  necessary  to  run  GIMAGE  several  times  on  a  given  sector  in  order  to 
properly  center  the  storm.)  The  resultant  image  is  then  processed  with  the 
program  TCIF  (for  Tropical  Cyclone  Intensity  Forecast)  to  better  estimate  the 
storm  center  location  (in  both  pixel/raster  and  latitude/longitude  coordinates). 
The  program  ZEBFLS  is  then  run  to  transfer  the  image  produced  by  GIMAGE  to 
magnetic  tape  in  a  format  acceptable  to  the  FNOC  computer  SPC.  Finally, 
ANGLEL  is  run  to  determine  the  amount  of  rotation  needed  for  each  image  in 
the  Lagrangian  transformation. 
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ORDER  OF  EVENTS 


ON  5PADS 


1.  Generate  large  sector  image  files  in  for¬ 
mat  acceptable  to  5PADS 


2.  5ubset  sector  image  files  to  generate 
512  x  512  files  centered  on  storm 


3.  Find  coordinates  of  spiral  centers  in  both 
(pixel,  raster)  and  (latitude,  longitude) 


4.  Store  512  x  512  images  on  tape  in  format 
acceptable  to  5PC 


5.  Calculate  rotation  angle  needed  for  each 
image 


ON  SPC 


6.  Rotate  each  image  and  save  new  images 
on  tape 

7.  Replace  missing  images  (if  any)  using 
interpolation  scheme  and  save  on  tape 


ON  SPADS 


8.  Read  new  images  from  tape  to  disk 

9.  Display  new  images  chronologically  and 
transfer  to  videotape  using  RBG  to 
NTSC  color  encoder 


FIGURE  2-1:  SUMMARY  OF  IMAGE  TRANSFORMATION  PROCESS 
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The  program  ANGLEL  and  its  subroutine  CKIN  are  listed  in  Appendix  A 
along  with  their  flow  charts.  This  program  calculates  the  angle  that  an  image 
must  be  rotated  in  order  to  put  it  into  the  Lagrangian  frame.  The  center  of  the 
storm  (determined  using  TCIF)  is  the  origin  of  the  Lagrangian  frame  and  the 
vertical  axis  becomes  the  direction  of  its  motion.  If  two  consecutive  images  are 
denoted  by  the  subscripts  i  and  i  +  1  then  the  change  in  latitude  and  longitude  of 
the  storm  center  is 


A  lat  =  lat.  .  -  lat. 

l+l  i 

A  long  =  l°ngi+1  -  long. 

The  angle  of  rotation,  0  ,  needed  to  put  the  image  into  the  Lagrangian  frame  is 

then  given  by 


tan  0  =  tan(cos(lat)  •  A  long! /sin(  A  lat) 

where  lat  is  the  higher  latitude  of  the  two  points.  In  order  to  speed 
computations,  this  equation  was  approximated  to 

tan  0  =  cos  (lat)  (A  long/A lat) 

since  A  lat  and  A  long  are  such  small  angles.  (The  inaccuracy  introduced  by  this 
approximation  is  considerably  smaller  than  the  inaccuracies  of  the  A  lat  and 
A  long  measurements  caused  by  errors  in  storm  center  positioning.)  If  the  second 
point  (i+1)  is  equatorward  of  the  first  point  (i)  then  0  is  adjusted  by  adding  or 
subtracting  180°.  The  results  of  using  this  program  are  discussed  in  section  3.0. 

2.2  Processing  on  SPC 

The  images  stored  on  tape  by  ZEBFLS  can  now  be  processed  with  the  SPC 
computer  at  FNOC.  As  indicated  in  Figure  2-1,  two  programs  are  used  in  this 
processing.  The  first  is  TCROTAT,  which  reads  an  image  from  tape,  rotates  it 
through  the  angle  calculated  by  ANGLEL,  and  writes  the  new  image  to  another 
tape.  The  second  is  ITRPIMG,  which  creates  an  interpolated  image  from  two 
other  images  (the  "before”  and  "after"  images  that  had  already  been  rotated) 
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and  writes  the  resultant  image  to  tape.  These  programs  are  executed  with  the 
job  card  streams  3TCR0TATE  and  3INTERPOLATE.  These  programs  will  now  be 
discussed  in  detail. 

2.2.1  Program  TCROTAT 

The  job  stream  listing,  flow  charts  and  program  and  subroutine  listings  of 
TCROTAT  are  given  in  Appendix  B.  The  input  data  to  this  program  consists  of 
the  storm  center  location  in  (pixel,  raster)  coordinates  (as  found  by  running 
TCIF)  and  the  angle  the  image  needs  to  be  rotated  (as  found  by  running 
ANGLEL).  Columns  1-3  of  the  input  record  contain  ICNTR,  the  pixel  value  of 
the  storm  center;  columns  5-7  contain  3CNTR,  the  raster  value  of  the  center; 
and  columns  9-15  contain  ANGLE,  the  rotation  angle  for  the  image  (positive 
counterclockwise).  (Warning:  Do  not  invert  pixel  and  raster  number  when  using 
TCROTAT.  This  is  easy  to  do  since  in  TCROTAT  the  pixel  is  labeled  I  and  the 
raster  3  whereas  in  TCIF  the  notation  is  reversed.)  The  user  also  has  the  option 
of  using  the  line  printer  to  display  the  central  portion  (128  x  128)  of  the  storm 
for  both  the  original  and  rotated  images;  this  option  is  exercised  through  the 
setting  of  switches  in  the  job  control  stream.  The  program  thus  starts  by  reading 
the  switches  and  the  input  data  (storm  center  and  rotation  angle)  and  checking 
for  validity. 

The  image  file  is  then  read  from  the  tape  (referred  to  as  TAPE4  in  the 
program)  and  stored  in  Extended  Core  Storage  (ECS).  An  image  file  consists  of 
515  records:  record  1  is  the  header,  records  2-513  are  image  data  (one  raster  per 
record),  and  records  514  and  515  are  documentation.  Each  raster  (or  record) 
contains  512  8-bit  pixels  (or  bytes)  packed  into  69  words;  the  image  file 
therefore  is  69  x  515  =  35535  words  long.  As  the  image  is  read  in,  it  is  shifted 
such  that  the  spiral  center  location  will  be  at  pixel  256  and  raster  256.  This 
means  that  blank  rows  and  columns  will  be  inserted  where  the  image  has  moved 
away  from  the  edges.  The  first  word  address  (FWA)  of  each  raster  is  stored  in 
IECl(l)  through  IECl(512),the  FWA  of  the  header  record  is  stored  in  IEC1(513) 
and  the  FWA's  of  the  documentation  are  stored  in  IEC1(514)  and  IEC1(515).  This 
array  is  contained  in  the  COMMON  block  EQBLK,  along  with  a  similar  array  for 
the  rotated  image.  This  second  array,  IEC2,  is  512  words  long,  having  only  the 
FWA's  of  the  rasters  of  the  rotated  image  since  the  header  and  documentation 
records  will  be  the  same  as  for  the  original  image. 
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If  switch  1  is  set,  the  central  portion  (128  x  128)  of  the  original  image 
(shifted  to  the  storm  center)  will  be  displayed,  using  the  line  printer  as  the 
display  devise.  Since  the  line  printer  outputs  characters  at  10  per  inch 
horizontally  and  6  or  8  per  inch  vertically,  the  resultant  image  is  elongated 
vertically.  This  distortion  is  of  no  consequence  since  these  intermediate  displays 
are  only  for  verification  of  storm  center  and  rotation. 

The  line  printer  display  is  executed  with  a  call  to  the  subroutine  IMGPLOT, 
which  is  passed  the  argument  ITYP.  If  ITYP  =  1,  the  original  image  is  displayed; 
if  ITYP  =  2,  the  rotated  image  is  displayed  on  the  line  printer.  This  subroutine 
reads  the  middle  128  rasters  of  the  image  from  ECS.  It  unpacks  the  data,  assigns 
a  printer  character  for  each  pixel  intensity  according  to  Table  2-1,  and  prints  the 
line.  As  stated  above,  only  the  central  128  x  128  portion  of  the  image  is 
displayed. 

The  image  is  then  rotated  with  a  call  to  subroutine  ROTATE,  which  is 
passed  the  argument  ANGLE.  The  latter  is  the  angle,  in  degrees,  the  image  is  to 
be  rotated,  with  positive  angles  being  counterclockwise.  In  order  to  save 
computer  time  the  corners  of  the  512  x  512  image  are  not  rotated,  since  they 
would  be  rotated  out  of  the  image  and  be  lost  anyway.  This  is  accomplished  by 
calculating  the  start  and  end  points  for  each  raster 


START  PT  =  512  -  END  PT 


where 


Y  =  raster  number  -  1 


If  (x,y)  is  the  location  of  a  pixel  in  the  original  image,  then  the  location 
(x',y')  of  that  pixel  in  the  rotated  image  can  be  calculated  from 


x'  =  x  •  cos(ANGLE)  +  y  •  sin(ANGLE) 
y'  =  -x  •  sin(ANGLE)  +  y  •  cos(ANGLE) 


If  I  is  the  pixel  number  (I  =  1  at  the  left  edge  of  the  screen)  then  x  =  I  -  256. 
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TABLE  2-1 


CORRESPONDENCE  BETWEEN  SATELLITE  IMAGE  INTENSITY 
RANGES  AND  LINE  PRINTER  DISPLAY  CHARACTERS 


INTENSITY  RANGE 

CHARACTER 

0-15 

BLANK 

16-31 

.  (PERIOD) 

32-47 

:  (COLON) 

48-63 

-  (HYPHEN) 

64-79 

+ 

80-95 

* 

96-111 

1 

112-127 

2 

128-143 

3 

144-159 

4 

160-175 

5 

176-191 

6 

192-207 

7 

208-223 

8 

224-239 

9 

240-255 

= 
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Likewise,  if  3  is  the  raster  number  (3  =  1  at  the  top  of  the  screen),  then  y  =  256  - 
J.  The  center  of  the  screen  is  therefore  at  (x,  y)  =  (o,  o)  and  (I,  J)  =  (256,  256). 

The  coordinates  (x,  y)  of  all  the  pixels  in  the  original  image  are  integers 

but  the  coordinates  (x',  y')  of  the  pixels  in  the  rotated  image  are,  in  general,  not 

integers.  A  bilinear  interpolation  was  therefore  performed  to  assign  the  proper 

intensities  to  the  integer  pixel  location  of  the  new  image.  If  x'  of  the  new  image 

falls  between  xi  and  x.+1  of  the  original  image,  and  if  y'  falls  between  y.  and 

v.  t , then 
3  j+r 

I(x',  y')  =  Io  +  (Ij  -  I0)  Ax  +  (I2  -  IQ)  Ay  +  (I3  +  IQ  -  I2  -  IL)  Ax  Ay 

where 

'o  ■ 1  (xi’  Vj> 

h  =  I<xi+i’V 
*2  =  1  <Xi’  V> 
b =  1  k'i+i’  yj+P 

Ax  =  x'  -  x. 

1 

Ay  =  y'  -  y- 

In  practice,  the  (x',  y')  are  kept  integer  and  the  corresponding  (x,  y)  locations 
calculated,  then  the  above  interpolation  equation  is  used  to  assign  the  proper 
intensity  to  the  pixel  at  (x',  y'). 

In  order  to  minimize  central  memory  (CM)  requirements,  the  new  image  is 
built  one  raster  at  a  time,  bringing  into  CM  only  those  portions  of  the  original 
image  that  are  necessary.  The  new  image  is  packed,  in  ECS,  in  the  same  format 
as  the  original  image. 

Upon  completion  of  image  rotation,  switch  2  is  tested.  If  it  is  ON  then 
IMGPLOT  is  called  again,  this  time  to  display  the  central  portion  of  the  rotated 
image. 
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Finally,  the  rotated  image  is  written  to  TAPE3  along  with  the  same  header 
and  documentation  as  on  the  original  image.  The  output  tape  has  a  format 
identical  to  that  of  input  tape  as  it  will  be  read  by  SPADS. 

In  using  3TCROTATE,  the  user  must  take  care  to  change  the  number  of 
files  skipped  on  each  tape  as  well  as  use  the  proper  VSN  (volume  serial  number) 
for  each  tape,  each  time  this  job  stream  is  executed. 


2.2.2  Program  ITRPIMG 

The  job  stream  listings,  flow  charts,  and  program  and  subroutine  listings  of 
ITRPIMG  are  given  in  Appendix  C.  This  program  replaces  missing  or  garbled 
images  by  interpolating  between  "before"  and  "after"  images.  Since  the  "before" 
and  "after"  images  generally  have  different  orientations  they  must  be  rotated 
slightly  relative  to  one  another  before  the  interpolation  is  done.  The  amount  of 
rotation  needed  is  determined  by  a  cross-correlation  calculation,  in  which  the 
pixel  intensities  on  a  circle  of  arbitrary  radius  (input  by  the  user)  on  one  image 
are  multiplied  by  the  corresponding  pixel  intensities  of  the  other  image.  The 
sum  of  these  products  is  a  measure  of  the  degree  of  cross-correlation  between 
the  two  images.  By  rotating  these  circles  relative  to  each  other  and  repeating 
the  cross-correlation  calculation  for  each  relative  orientation,  the  "best"  orien¬ 
tation  is  determined  by  noting  where  the  cross  correlation  is  a  maximum.  This 
process  is  shown  schematically  in  Figure  2-2. 

The  program  starts  by  reading  switch  1  and  the  input  "card",  which 
contains  the  user-supplied  radius  (13  format)  of  the  circle  used  for  the  cross¬ 
correlation  calculation.  This  radius,  which  is  measured  in  the  unit  corresponding 
to  the  size  of  the  pixel,  must  be  positive  and  less  than  255.  The  program  verifies 
this,  then  reads  in  the  "before"  and  "after"  images  with  a  call  to  the  subroutine 
READIMG;  the  "before"  image  is  read  from  TAPE3  and  the  "after"  from  TAPE4. 
(In  the  job  stream  3INTERPOLATE  be  sure  the  correct  VSN  is  specified  and  that 
the  proper  files  have  been  copied  to  TAPE3  and  TAPE4  prior  to  executing 
ITRPRMG.) 
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Start  with  "before"  and  "after" 
images.  User  specifies  radius  of 
circle  for  cross-correlation. 


Determine  intensities  around  cir¬ 
cumference  of  each  cross-corre¬ 
lation  circle,  interpolating  where 
necessary.  The  step  size,  A  9, 
corresponds  to  one  pixel  dia¬ 
meter. 


J 

Image  A 


Image  B 


Calculate  cross-correlation  bet¬ 
ween  1^  (  9)  and  Ig  (  0) 

VIV'V'b'V"  as) 

i 

for  various  values  of  n  A  0  and 
pick  value  of  n  giving  highest  S  n  . 


Calculate  interpolated  image  by 
rotating  image  A  by  Yi  nmax  A  0 
and  image  B  by  -h  nmax  A  9  and 
averaging  each  pixel. 


Interpolated 

Image 


FIGURE  2-2:  INTERPOLATION  PROCESS 


The  images  are  stored  in  ECS.  The  first  word  addresses  (FWA)  are  stored  in 
IEC1  and  IEC2.  IEC1  being  the  FWA's  of  the  before  image  and  IEC2  being  the 
FWA's  of  the  after  image.  The  images  are  stored  in  the  same  format  as  in 
TCROTAT. 

If  switch  1  is  ON,  the  central  portions  (128  x  128  pixels)  of  each  image 
("before"  and  "after")  are  displayed,  using  the  line  printer,  by  a  call  to  the 
subroutine  PRT.  This  subroutine  is  almost  identical  to  the  subroutine  IMGPLOT, 
which  is  associated  with  the  program  TCROTAT. 

The  program  then  determines  the  number  of  pixels  that  would  fit  on  the 
cross-correlation  circle  (this  equals  the  circumference)  and  calculates  the 
angular  displacement,  A0  >  that  corresponds  to  the  width  of  one  pixel  at  the 
distance  r,  the  radius  of  the  circle.  The  (x,y)  coordinates  of  each  pixel  on  the 
circle  in  steps  of  A0  are  then  determined  and  the  corresponding  pixel  intensities 
found  with  a  bilinear  interpolation;  this  is  done  for  both  images.  The  cross¬ 
correlation  of  the  resulting  one-dimensional  intensity  arrays  is  then  given  by 

V  I  >A  (©i’  fe  (Si+n  AS) 

i 

where  the  subscripts  A  and  B  refer  to  the  two  images,  i  is  the  pixel  number,  and 

n  is  the  measure  of  the  amount  one  array  is  shifted  relative  to  the  other.  The 

maximum  correlation  corresponds  to  the  largest  S  .  To  save  computer  time,  the 

S  are  first  calculated  from  n  =  -25  to  n  =  +25  in  steps  of  a  n=5.  Then  n  ,  the 

value  of  n  for  the  maximum  Sn  from  these  coarse  resolution  calculations,  is  used 

as  the  mid-point  of  fine  resolution  (An=l)  calculations  of  S  over  the  interval 

n 

n  ±5.  Then  n  .  the  value  of  n  for  the  greatest  S  within  this  range,  is 
m  max  °  n  07 

determined.  The  needed  rotation  of  one  image  relative  to  the  other  to  achieve 
maximum  correlation  is  then  n 

max 

The  above  mentioned  bilinear  interpolation  is  performed  with  a  call  to 
subroutine  INTRP.  The  arguments  passed  to  this  subroutine  by  the  calling 
program  are  IF1LE,  the  number  of  the  image  file  (1  for  "before"  or  2  for  "after") 
and  the  pixel  and  raster  numbers  (x,y).  The  subroutine  then  reads  from  ECS  the 
raster  on  either  side  of  y  and  calculates  the  intensity  at  pixel  location  x  with  a 
bilinear  interpolation;  this  intensity  is  returned  to  the  calling  program  through 
the  argument  INT. 
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The  cross-correlation  calculation  and  determination  of  the  nmor  nmax  is 
performed  with  a  call  to  subroutine  SUMSN.  The  arguments  passed  to  this 
subroutine  are  NOPTS,  the  number  of  points  (pixels)  on  the  cross-correlation 
circle  and  NINC,  the  increment  in  n  to  use.  The  subroutine  returns  the 
argument  NMAX,  the  value  of  n  the  corresponds  to  the  maximum  Sn  for  this  NINC. 

The  program  then  rotates  each  image  by  half  of  the  total  required:  the 
first  image  is  rotated  Yi  nmax  A  0  and  the  second  nmax  A0-  The  rotated 
images  are  then  averaged,  pixel  by  pixel,  to  generate  the  final  image,  which  is 
written  to  ECS.  These  calculations  are  performed  through  a  call  to  subroutine 
ROTAVG,  which  is  passed  the  argument  ANGLE  (=  Yi  r»maxA  ©)• 

Switch  1  is  then  tested  again  and  if  found  to  be  ON,  the  central  portion  of 
the  new  image  is  displayed  with  the  line  printer  with  a  call  to  subroutine  PRT, 
described  earlier. 

Finally,  the  new  image  is  written  to  TAPE5  with  a  call  to  subroutine 
RITIMG.  The  header  record  is  first  output;  it  is  generated  as  "IMAGE 
INTERPOLATED  AFTER"  followed  by  the  first  15  characters  taken  from  the 
header  of  the  "before"  image.  The  rasters  are  then  read  from  ECS  one  at  a  time 
and  output  to  TAPE5.  Finally,  the  two  documentation  records  from  the  "before" 
image  are  copied  to  TAPE5.  This  file  is  written  to  magnetic  tape  in  the  job 
stream  following  execution  of  the  program.  (Be  sure  the  proper  VSN  is  used  and 
that  the  proper  number  of  files  have  been  skipped  prior  to  copying  TAPE5  to 
magnetic  tape.) 

2.3  Final  Processing  on  SPADS 

After  all  of  the  rotated  images  have  been  generated  and  saved  on  tape  at 
SPC,  they  are  ready  for  display  on  the  SPADS  and  transfer  to  videotape.  The 
images  are  transferred  to  videotape  in  chronological  order  with  a  one  second 
view  time  of  each  frame;  on  playback  this  gives  the  appearance  of  time-lapse 
photography.  Both  the  original  images  (which  are  unrotated  and  not  perfectly 
centered)  and  the  new  ones  (those  in  the  Lagrangian  frame)  have  been  saved  on 
videotape.  In  section  3  we  discuss  these  results. 
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The  hardware  needed  for  this  process  is  an  RBG  to  NTSC  color  encoder. 
The  Genisco  display  unit  outputs  digital  red-blue-green  (RBG)  signals  but  since 
TV  monitors  and  video  tape  recorders  need  to  receive  NTSC  type  signals,  there  is 
a  need  to  convert  the  signals  from  one  type  to  the  other.  ODSI  personnel  used 
the  LENCO  CCE-850  unit  for  this  project.  This  unit  is  ODSI  property,  so  the 
government  will  have  to  obtain  its  own  LENCO  unit  if  it  wishes  to  repeat  this 
process  on  another  cyclone,  or  follow  the  alternative  procedure  outlined  at  the 
end  of  this  section. 

A  schematic  diagram  for  the  attachment  of  the  LENCO  unit  to  the  Genisco 
image  display  processor  is  shown  in  Figure  2-3.  As  can  be  seen,  the  attachment 
is  quite  simple.  Not  shown  in  the  figure  is  the  optional  attachment  of  a  color  TV 
monitor  to  the  video  tape  recorder.  Even  though  the  use  of  the  TV  monitor  is 
not  absolutely  necessary,  it  is  advisable  to  use  one  during  the  recording  session 
so  that  the  color  adjustments  on  the  LENCO  Unit  can  be  optimized. 

One  additional  hardware  modification  is  needed  to  complete  this  process. 
Since  the  Genisco  image  display  controller  normally  outputs  image  frames  at 
the  rate  of  40  per  second,  and  the  TV  monitor  and  video  recorder  need  an  input 
rate  of  30  frames  per  second,  a  serious  synchronization  problem  arises.  To 
correct  this,  the  video  output  board  in  slot  3  of  the  PGP  unit  must  be  replaced 
with  one  that  outputs  the  needed  30  frames  per  second.  The  proper  board  was 
obtained  on  loan  from  Genisco  for  the  recording  session.  The  high  resolution 
RBG  color  monitor  becomes  useless  when  this  change  is  made,  as  it  needs  input 
at  a  rate  of  40  frames  per  second.  This  is  one  more  reason  for  using  a  TV 
monitor  during  the  recording  session. 

The  SPADS  program  5DISPLAY  was  used  for  the  recording  session;  any  of  a 
number  of  other  routines  could  have  been  used,  however.  The  procedure  is  to 
first  load  the  entire  sequence  of  images  from  tape  onto  disc.  For  the  original 

I 

images,  this  is  accomplished  with  the  LOAD  command  but  for  the  new  images 
the  program  ZEBFLS  must  be  used.  The  videotape  first  is  allowed  to  record 
several  seconds  of  blank  image.  The  tape  is  then  stopped  with  the  PAUSE 
button.  Images  are  then  brought  to  the  screen  in  chronological  sequence  and  a 
one-second  recording  of  each  image  is  made,  again  using  the  PAUSE  button.  It  is 
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important  to  use  the  PAUSE  button  instead  of  RECORD  and  STOP  since  in  the 
latter  case  the  tape  motor  has  to  come  up  to  speed  before  recording  can  start, 
which  causes  blanks  to  occur  in  the  sequence;  this  also  makes  it  difficult  to 
estimate  the  proper  length  of  time  to  record. 

Following  the  sequence  of  original  images,  several  seconds  of  blank  image 
were  recorded,  followed  by  the  sequence  of  new  images  in  the  Lagrangian  frame. 
In  using  5DISPLAY  we  chose  to  color  slice  in  the  temperature  intervals  0°,  -12°, 
-24°,  -36°,  -48°,  -60°  and  -7 2°C.  The  government  may  wish  to  repeat  the 
recording  with  another  temperature  set  in  order  to  emphasize  different  cloud 
features. 

An  alternative  procedure  to  using  the  LENCO  unit  is  to  photograph  each 
image  on  the  high  resolution  color  monitor.  Color  slides  can  then  be  copied  to 
videotape  using  the  services  of  a  commercial  photography  company  such  as 
FOTOMAT.  (Last  year  they  were  offering  this  service  at  a  rate  of  $9.00  per  100 
35mm  slides).  With  this  procedure,  care  must  be  taken  to  not  disturb  the  camera 
when  changing  film.  An  advantage  of  this  technique  is  that  a  label  with  the  day 
and  hour  (such  as  IGNACI0296.1 8)  could  be  affixed  to  the  screen  prior  to 
photographing  the  image,  thus  providing  a  "clock"  on  the  image.  Also,  titles 
could  be  inserted  to  introduce  the  sequence.  A  disadvantage  is  that  film 
processors  take  little  care  in  framing  slides  so  that  the  images  will  probably  not 
be  properly  centered.  With  custom  processing,  however,  this  may  not  be  a 
problem. 
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3.0  DISCUSSION  OF  RESULTS 


The  images  of  Hurricane  Ignacio,  as  captured  in  sequence  on  videotape 
in  both  the  original  form  and  in  Lagrangian  coordinates,  reveal  a  number  of 
interesting  meteorological  features  in  cloud  morphology  and  motions.  Study 
of  these  sequences  also  reveals  some  problems  that  should  be  corrected  in 
future  work;  these  problems  could  not  be  corrected  in  this  study  due  to  time 
limitations,  as  explained  below. 

3.1  Motions  in  the  Eulerian  Frame 


The  Government,  during  work  on  another  project,  used  program  GIMAGE 
to  subset  the  large  sector  image  files  of  Hurricane  Ignacio  and  produce  approximately 
centered  images  of  this  storm  that  were  saved  on  tape.  The  government  also 
ran  program  TCIF  (for  another  project)  to  determine  more  accurate  locations 
of  the  storm  centers.  These  images  and  centers  were  provided  to  ODSI  as 
the  basic  data  set  for  this  contract.  The  storm  track,  based  on  these  data, 
is  shown  in  Figure  3-1;  the  gaps  marked  "M"  denote  missing  or  garbled  images. 
These  images  have  been  transferred  to  videotape  in  chronological  sequence, 
as  described  above.  The  resulting  video  sequence  is  somewhat  "jumpy"  since 
these  images  are  not  perfectly  centered  (something  that  is  difficult  to  do  with 
GIMAGE),  but  nevertheless  the  sequence  shows  some  interesting  features. 

The  most  obvious  motion  that  can  be  seen  in  this  first  sequence  is  the 
counter-clockwise  rotation  of  the  central  portion  of  the  storm  throughout  its 
life.  The  central  dense  overcast  (CDO)  portion  of  the  storm  consists  of  convective 
clouds  with  great  vertical  development.  Although  motions  within  the  CDO 
cannot  be  discerned  with  the  temperature  slicing  used  in  this  sequence,  the 
motions  along  the  edges  of  the  CDO  are  obvious.  Large  convective  cells  are 
seen  to  form  and  dissolve  at  a  distance  of  100  to  200  miles  from  the  storm 
center.  These  cells,  which  are  sometimes  larger  than  the  CDO,  move  counter¬ 
clockwise  about  the  CDO  very  slowly;  at  times  they  seem  to  compete  with 
the  CDO  for  becoming  the  new  center  of  the  cyclone.  The  outer  spiral  banded 
structure  (OSBS)  sometimes  rotates  about  the  CDO  and  at  other  times  remains 
stationary. 
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HURRICANE  IGNACIO,  1979 
STORM  TRACK 


LONGITUDE 


FIGURE  3-1:  THE  STORM  TRACK  OF  HURRICANE  IGNACIO 
The  gaps  marked  M  indicate  missing  or  garbled  images.  The  X's  mark  the 
position  of  the  storm  at  the  beginning  of  each  day;  the  corresponding  numbers 
indicate  the  day  and  hour,  where  the  day  number  is  counted  from  January  1, 
1979. 
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The  size  of  the  storm  varied  considerably  during  its  life.  At  the  start 
of  the  sequence  the  storm  fills  the  entire  screen;  that  is,  it  had  a  diameter 
of  over  1000  miles;  the  CDO  had  a  diameter  of  perhaps  one-tenth  this  size. 

In  the  middle  of  the  sequence  the  cyclone  had  shrunk  to  half  of  its  original 
size  and  had  a  well-developed  "bullseye"  type  central  eye.  Near  the  end  of 
the  sequence  the  storm  had  grown  again  and  nearly  attained  its  original  size. 
During  the  period  of  its  smaller  size,  the  OSBS  was  very  sharply  defined  but 

at  the  other  times,  it  was  more  patchy. 

3.2  Motion  in  the  Lagrangian  Frame 

Program  ANGLEL  was  used  to  calculate  the  rotation  angles  needed  to 
produce  images  in  the  Lagrangian  frame.  The  results  are  given  in  Table  3- 
1  and  shown  in  Figure  3-2.  These  angles  were  used  as  input  data  for  program 
TCROTAT  for  producing  the  second  sequence  of  images  seen  on  the  videotape. 

Samples  of  the  images  generated  by  TCROTAT  are  shown  in  Figures  3-3 
and  3-5.  A  sample  of  the  image  interpolated  between  those  two  images  is  shown 
in  Figure  3-4;  the  latter  was  generated  with  program  ITRPIMG.  It  was  originally 
thought  that  the  image  of  Ignacio  for  day  298,  hour  12,  (IGNACIO  298.12)  was 
unavailable  due  to  a  bad  sector  scan  number  (that  is  why  the  interpolated  image 
shown  in  Figure  3-4  was  generated),  but  it  was  later  discovered  that  IGNACIO 
298.12  was  readable;  it  is  showm  in  Figure  3-6.  The  image  in  the  figure  has  been 
rotated  by  hand  when  photocopying  so  that  its  orientation  is  the  same  as  that  of 
the  interpolated  image.  Comparison  of  Figures  3-4  and  3-6  shows  that  image 
generated  by  ITRPIMG  is  very  good  approximation  to  the  true  image,  the 
corresponding  cross-corrrelation  function  is  shown  in  Figure  3-7;  its  peak  is 
broad  but  smooth,  and  the  maximum  is  well  defined. 


These  Lagrangian  images  show  features  not  obviously  in  the  first  sequence. 
Perhaps  the  most  notable  is  the  tendency  for  the  large,  well  developed  convective 
cells  to  lie  in  the  direction  of  motion  of  the  storm,  or  somewhat  to  the  right 
of  this  direction;  another  aspect  of  this  observation  is  that  these  convective 
cells  do  not  seem  to  rotate  about  the  CDO,  as  they  did  in  the  first  sequence. 
Early  in  the  life  of  the  storm  rotation  is  absent  in  the  Lagrangian  frame.  Later, 
however,  a  pronounced  rotation  does  occur.  Late  in  the  storm  a  large  convective 
cell  forms  in  the  OSBS;  significantly,  this  cell  also  lies  ahead  of  the  CDO, 
that  is,  in  the  direction  of  motion  of  the  storm. 
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TABLE  3-1:  LAGRANGIAN  ROTATION  ANGLES 


TIME 

ANGLE 

DAY  HR 

TO  ROTATE 

296  .  1  8 

-112.71 

296.20 

-  109 . 42 

296 . 22 

-110.63 

297.00 

-  105 . 24 

297.02 

-107.88 

297.04 

-111.26 

297.06 

-97 . 27 

297.08 

-86.56 

297  .  10 

-  103 . 66 

297  .  1  2 

-85  .  14 

297 . 14 

-73 . 28 

297  .  1  6 

-66.94 

297 . 20 

-62 . 93 

297.22 

-90.00 

298.01 

-90.00 

298.02 

-92  .  25 

298 .04 

-77.65 

298.06 

-66 . 34 

298 .08 

-69.02 

298  .  10 

-59.04 

298.14 

-70.61 

298  .  1  6 

-52  .  17 

298  .  1  8 

-64.43 

298 . 20 

-62 . 84 

298 . 22 

-60.76 

299 .00 

-55 . 60 

299.02 

-28 . 92 

299 .04 

-28.21 

299 . 06 

-39.71 

299.08 

-31.35 

299.10 

-55 .79 

299 . 12 

-46.60 

299  .  1  4 

-52 . 54 

299.  1  6 

-38.94 

299 . 1 8 

-61.63 

299.20 

-49.10 

299.22 

-65.  19 

300.00 

-72 . 66 

300.02 

-51.17 

300.04 

-56 . 72 

300.06 

-39 . 80 

300.08 

-61.42 

300. 10 

-92 . 39 

300. 12 

-54  .  65 

300. 1 4 

-64  .  23 

300. 1 6 

-120.65 

300. 1 8 

-126.73 

300.20 

-78 .21 

300.22 

-43 . 73 

301 .01 

-12.44 

301 .02 

-35 .01 

301.04 

-120.18 

301 .06 

-20.37 

301 .08 

26.79 

301 . 10 

-7.74 

301.12 

.00 

301.14 

-130.06 

301 . 1 6 

81  .05 

301.20 

82 . 52 

301.22 

136.41 

302.00 

-96 . 66 

302.02 

43 . 59 

302 .04 

10.76 

302.06 

88.49 

302.08 

83 . 33 

302 . 10 

111.85 

302. 12 

108.38 

302. 14 

105.19 

302. 16 

77 . 52 

302 . 1 8 

32 .32 

302.20 

87 .69 

302 .22 

92.08 

303.00 

108.37 

303.02 

96  .  97 

303.04 

102 . 26 

303.06 

102 . 26 

LATITUDE  LONGITUDE 


12. 

15 

-96  . 

24 

1  1  . 

97 

-96  . 

68 

1  1  . 

87 

-96  . 

97 

1  1  . 

80 

-97  . 

1  6 

1  1  . 

76 

-97  . 

3  1 

1  1  . 

70 

-97  . 

50 

1  1  . 

62 

-97  . 

71 

1  1  . 

59 

-97  . 

95 

1  1  . 

60 

-98  . 

12 

1  1  . 

55 

-98  . 

33 

1  1  . 

57 

-98  . 

57 

1  1  . 

67 

-98  . 

9  1 

1  1  . 

82 

-99  . 

27 

1  1  . 

92 

-99  . 

47 

1  1  . 

92 

-99  . 

60 

1  1  . 

92 

-99  . 

8  1 

1  1  . 

9  1 

-100. 

07 

1  1  . 

94 

-100. 

21 

12  . 

00 

-100. 

35 

12  . 

1  2 

-  100. 

67 

12  . 

29 

-  100. 

96 

12. 

40 

-101  . 

28 

12. 

65 

-101  . 

61 

12  . 

79 

-  101  . 

9  1 

12. 

92 

-  102  . 

17 

1  3  . 

04 

-102  . 

39 

1  3  . 

1  8 

-  102  . 

60 

1  3  . 

69 

-  102  . 

89 

1  4  . 

1  6 

-103  . 

15 

1  4  . 

37 

-  103  . 

33 

1  4  . 

64 

-  103  . 

50 

1  4  . 

87 

-  103  . 

85 

1  5  . 

08 

-  104  . 

08 

15  . 

25 

-  104  . 

31 

1  5  . 

56 

-104  . 

,  57 

15  . 

,69 

-  104  . 

,  82 

1  5  . 

,  84 

-105  . 

.00 

1  5  . 

.  96 

-105  . 

,  27 

1  6  . 

.05 

-  105  . 

.  57 

1  6  , 

.  22 

-105  . 

.79 

16  , 

.  39 

-  106  , 

.06 

1  6 

.62 

-106  , 

.  26 

1  6 

.74 

-  106  , 

.  49 

1  6 

.73 

-106 

.  74 

16 

.  92 

-107 

.02 

1  7 

.04 

-  107 

.  28 

1  6 

.87 

-107 

.  58 

1  6 

.77 

-107 

.  72 

1  6 

.  78 

-  107 

.  77 

1  6 

,91 

-107 

.  90 

17 

.04 

-107 

.93 

17 

.  1  9 

-  108 

.04 

17 

.  1  4 

-  108 

.  13 

17 

.  32 

-108 

.  20 

17 

.  49 

-  108 

.  1  1 

17 

.  84 

-  108 

.  1  6 

17 

.95 

-108 

.  1  6 

17 

.  83 

-108 

.31 

17 

.  86 

-108 

.  1  1 

17 

.  87 

-108 

.03 

17 

.  77 

-  107 

.  93 

17 

.76 

-  108 

.02 

17 

.  82 

-  107 

.  96 

18 

.  1  7 

-  107 

.  89 

18 

.  1  8 

-107 

.  49 

18 

.21 

-107 

.  22 

1  8 

.  13 

-107 

.01 

18 

.07 

-  106 

.82 

17 

.  99 

••-106 

.51 

1  8 

.03 

-106 

.  32 

1  8 

.08 

-  105 

.  93 

1  8 

.09 

-105 

.  67 

1  8 

.08 

-  105 

.  38 

1  7 

.96 

-105 

.00 

1  7 

.  91 

-  104 

.57 

17 

.  85 

-104 

.  28 
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ROTATION  ANGLE  (DEGREES) 


HURRICANE  IGNACIO,  1979 


FIGURE  3-2:  LAGRANGIAN  ROTATION  ANGLES  FOR  HURRICANE  IGNACIO 
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NEW298. 10 


FIGURE  3-3 


The  image  produced  by  TCROTAT  from  IGNACIO  298.10.  From  Table  3-1 
we  see  that  the  rotation  angle  for  this  image  was  -59. 04°’ 


3-6 


FIGURE  3-4 


The  image  produced  by  ITRPIMG  by  interpolating  between  the  "before"  and 
"after"  images  shown  in  Figures  3-3  and  3-5. 
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The  image  produced  by  TCROTAT  from  IGNACIO  298.14.  From  Table  3-1 
we  see  that  the  rotation  angle  for  this  image  was  -70.61°. 
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FIGURE  3-6 


Image  IGNACIO 
to  give  the  same 


298.12;  the  original  image  has  been  rotated  during  photocopy 
orientation  as  Figure  3-4,  with  which  it  should  be  compared. 


Cross  correlation  of 
Ignacio  298.10  and  Ignacio  298.14 


Sn 


29545000  - 

29535000  - 

29525000  - 

29515000  “ 

29505000  - 

-42 


-36  -34  -32  -30 


N 


FIGURE  3-7:  THE  CROSS -CORRELATION  FUNCTION,  S  , 

AS  A  FUNCTION  OF  THE  NUMBER  OF  ANGULAR  n 

STEPS,  N,  FOR  THE  IMAGES  SHOWN  IN  FIGURES  3-3 

AND  3-5. 

The  radius  chosen  for  this  calculation  was  2 00  pixels.  The  peak 
occurs  at  N=-36,  which  corresponds  to  an  angle  of  -10.32°;  each 
image  is  thus  rotated  half  this  amount,  or  -5.16°,  prior  to  pixel- 
by-pixel  interpelation. 


-26 
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3.3  Recommendations  for  Future  Work 


As  interesting  as  the  two  sequences  on  videotape  are,  they  do  have  a 
number  of  defects  that  should  be  corrected  in  future  work.  Below  we  list  the 
most  obvious  shortcomings  and  suggest  ways  to  overcome  them.  These  defects 
were  not  corrected  in  the  present  contract  because  of  time  limitations.  The 
biggest  time  problem  is  that  of  processing  the  images  on  SPC;  usually  not 
more  than  two  or  three  images  could  be  rotated  in  one  day. 

3.3.1  Center  Locating 

Errors  in  the  placement  of  the  storm  centers  can  lead  to  serious  errors 
in  the  calculated  rotation  angles  for  the  images,  especially  when  the  storm 
is  moving  slowly.  This  can  be  seen  in  Figure  3-2;  during  days  300  and  301, 
especially,  the  calculated  angles  showed  a  lot  of  scatter.  Consequently,  many 
of  the  images  in  the  second  sequence  on  videotape  are  incorrect  and  actually 
confuse,  rather  than  help,  the  study  of  the  hurricane. 


Another  problem  caused  by  inaccurate  centering  is  in  the  interpolation 
process  to  replace  missing  images.  The  cross-correlation  technique  is  based 
on  the  assumption  that  consecutive  images  in  the  Lagrangian  frame  will  differ 
more  in  orientation  than  in  overall  structure.  If  the  centers  are  misplaced 
on  these  images,  however,  then  the  cross-correleation  will  probably  not  give 
an  accurate  value  of  the  relative  orientation;  hence,  the  interpolated  image 
will  not  be  as  good  as  it  should  be. 

To  determine  accurate  storm  centers  we  recommend  the  following  interative 
procedure. 


1.  Use  program  TCIF  to  establish  the  first-guess  centers  as  was  done 
for  this  project. 

2.  Calculate  both  the  velocity  of  the  storm  and  the  Lagrangian  rotation 
angle  for  each  image.  This  would  require  modifying  program  ANGLEL 
to  include  the  velocity  calculation. 
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3.  Plot  time  sequences  of  both  angle  and  velocity  and  then  smooth. 
Hand  smoothing  would  probably  be  accurate  enough  for  these  purposes, 
but  if  not  a  weighted  5  point  or  7  point  running  average  could  be 
used. 

4.  Use  the  smoothed  values  of  angle  and  velocity  to  recompute  the 
storm  track,  then  calculating  new  values  of  latitude/longitude  of 
the  storm  center  for  each  image.  Program  TCIF  could  then  be  run 
again  to  see  if  these  values  are  reasonable;  if  not,  adjustments  can 
be  made  to  establish  a  better  center  location. 

5.  Repeat  steps  2  through  4  as  many  times  as  necessary  to  obtain  a  self- 
consistent  set  of  angles  and  velocities  that  are  acceptable  to  the  user 
of  program  TCIF.  Special  attention  should  be  paid  to  slow  movement 
periods  and  to  turning  periods.  A  consistent  rotation  angle  variation 
must  be  followed. 

With  the  rotation  angles  obtained  this  way  the  Lagrangian  sequence  of 
images  would  be  more  meaningful. 

3.3.2  Additional  Screen  Information 


When  viewing  the  image  sequences  on  videotape,  it  would  be  very  useful 
to  have  additional  information  on  the  screen.  This  would  include: 

1.  The  day  and  hour  of  the  image; 

2.  The  velocity  of  the  storm; 

3.  The  Lagrangian  rotation  angle; 

4.  The  intensity  (maximum  wind  velocity)  of  the  cyclone,  as  obtained 
from  other  sources. 

It  would  also  be  useful  to  run  these  sequences  many  times  with  different 
temperature  slicing  in  order  to  emphasize  different  types  of  clouds  and,  in  at 
least  one  sequence,  to  show  the  temperature  structure  of  the  underlying  sea 
surface,  where  it  is  visible  between  the  clouds. 
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This  would  require  the  development  of  new  interactive  software  for 
the  input  of  labels  and  the  setting  of  video  lockup  tables. 

3.3.3  Mercator  Projection 


Hurricane  Ignacio,  as  viewed  from  the  GOES-West  satellite,  suffered  some 
distortion  in  its  images  since  it  was  located  about  30°  in  longitude  and  15°  in 
latitude  from  the  subsatellite  point.  Other  hurricanes  may  suffer  even  more 
distortion.  It  would  be  useful  to  process  the  large  sector  image  files  with  the 
program  GOESMRC  prior  to  processing  with  the  programs  discussed  in  this 
report  so  that  geometrical  distortions  could  be  removed  as  well  as  possible. 
Fortunately,  Tropical  Cyclones  are  confined  for  the  most  part  to  tropical 
latitudes  and  so  the  distortion  inherent  in  the  mercator  projection  are  not 

serious. 


3-13 


APPENDIX  A 


1.  Flow  chart  for  Program  ANGLEL 

2.  Listing  of  ANGLEL 

3.  Flow  chart  for  Subroutine  CKIN 

4.  Listing  of  CKIN 


A-l 


FLOW  CHART 
FOR  PROGRAM  ANGLEL 


A-2 


oooooononnoon 


THIS  PROGRAM  CALCULATES  THE  ANGLE  TO  ROTATE 
CYCLONE  IMAGE  ACCORDING  TO  THE  DIRECTION  OF 
MOVING  . 


A  TROPICAL 
WHICH  IT  IS 


INPUT: 

XLAT , XLON 

TIME 


-  MAX.  160  CENTER  LOCATIONS  OF  CYCLONE 

(-9999.  AFTER  LAST  POSITION) 

-  TIME  OF  POSITION  OF  CYCLONE 


OUTPUT 


THETA  -  MAX  160  ANGLES 

TIME  -  TIME  OF  POSITION 


OF  T. CYCLONE 


DIMENSION  XLAT ( 160) 
LOGICAL  OK 


XLON ( 1 60 ) ,  THETA( 160),  T I  ME ( 160) 


•PROGRAM  ANGLEL"  ,  „  „„ 

"THIS  PROGRAM  CALCULATES  THE  ANGLE  TO  ROTOTE 
“A  TROPICAL  STORM."  . 

"ENTER  TIME  OF  STORM  (DAY. HR  I.E.  299.02), 
"LATITUDE  AND  LONGITUDE  POSITION  AT  THE  PROMPT, 
"WHEN  NO  MORE  DATA  ENTER  -9999.  AT  THE  PROMPTS, 
=  0 


TYPE 
TYPE 
TYPE 
TYPE 
TYPE 
TYPE 
LAST 
0  =  1 

RADIAN  =  0.017453 

D°  10  1  ACCEPT  "ENTER  TIME:  "  , T I  ME ( 0  )  ,  " E NTER  L AT  I TUD E : " , X L AT ( J ) 
"ENTER  LONGITIUDE: "  ,XLON(0  )  , "<NL>" 

IF  (XLAT(O) .EQ.-9999. )  GOTO  20 

CALL  CKIN(TIME(U ) ,XLAT(0 ) , X L ON ( 0 ) , OK ) 

IF  (.NOT. OK)  GOTO  10 
0  =  0  +  1 

LAST  =  LAST  +  1 


10 
C 

C  CALCULATE 
C 


CONTINUE 

THETA  USING 


SLOPE  BETWEEN  POINTS 


20 


22 


25 

32 

30 


PRINT  2 
IF  (LAST 


IF 

IF 

DO 


EQ.160  .AND.  XLAT(LAST)  .NE.-9999, 
(LAST. LE. 160  .AND.  XLAT ( L AST  )  . EO . - 9 9 9 9 . 
(LAST. LE. 160  .AND.  XLAT( LAST ). NE . -9999 
30  L=  1  ,  LL 

DLAT  -  XLAT ( L  +  1 )  - 
DLON  =  XLON(L+l)  - 
I F ( DLAT . NE .0. ) 

IF ( DLON.GT.0. ) 

I F ( DLON . LT .0. ) 

GO  TO  32 
HL  -  XLAT ( L  +  1 > 

IF  (XLAT(L ) .GT.XLAT (L  +  l ) )  HL  =XLAT ( L  > 

A  «  COS ( HL*RAD IAN )  *  DLON/DLAT 

THETARAD  =  ATAN(A) 

THETA ( L  )  -  THETA RAD/ RAD  IAN 

IF  (XLAT(L+1  )  .GT. XLAT ( L  )  )  GOTO  32 
IF (THETA(L  )  . GT . 0 .  )  GO  TO  25 


L  L  • 
LL- 
LL- 


159 

LAST-2 

-LAST-1 


XLAT ( L ) 

XLON ( L ) 

GO  TO  22 
THETA ( L  >  =  90. 
THETA ( L )  =  -90. 


STORM  MOVE  TOWARD  EQUATOR 


THETA( L  ) 
GO  TO  32 
THETA ( L  ) 
PRINT  3 . 
CONTINUE 


THETA(L)  +  180. 

-  THETA ( L )  -  180. 
TIME(L),  THETA ( L ) 


XLAT ( L ) ,  XLON(L) 


C 

C 

2 


LL  =  LL  +  1 

THETA ( L L )  =  THETA <  L L -  1  > 

PRINT  3,  TIME(LL),  THETA(LL),  XLAT(LL),  XLON ( LL  ) 

GOTO  100 

FORMATS 

FORMAT ( " 1  TIME" ,T14 , "ANGLE" ,T22, "LATITUDE" ,T31 , "LONGITUDE"/ 
S  "  DAY  HR  " , T 1 2 , "TO  ROTATE"/) 

FORMAT! 1X.F7.2.3F10.2) 


3 
C 

C  ERROR  EXIT 
C 

100  STOP 

END 
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FLOW  CHART 
FOR  SUBROUTINE  CKIN 
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o  o  o  n  o  o  o 


SUBROUTINE  CK I N ( TIME , XLAT  * XLON , OK , F IRST ) 

THIS  SUBROUTINE  CHECKS  THE  VALIDITY  OF  THE  INPUT  DATA 

INPUT: 

TIME  -  TIME  OF  POSITION  OF  STORM 

XLAT,  XLON  -  LATITUDE,  LONGITUDE  POSITION  OF  STORM 


C 


10 


20 


30 


40 


STATIC  LASTME/0 .0/ 

INTEGER  TIME 
LOGICAL  OK 

OK  =  .TRUE. 

GOTO  20 

IF  (LASTME.LT. TIME)  GOTO  20 
OK  =  .FALSE. 

TYPE  “TIME  INVALID,  MUST  BE  LATER  THAN  PREVIOUS 
IF  (XLAT.GE.0.  .AND.  XLAT. LE. 25.)  GOTO  30 
OK  =  .FALSE. 

TYPE  “INVALID  LATITUDE  MUST  BE  0  TO  25., TRY  AGA 
IF  ( XLON . LE . -90.  .AND.  XL ON . G E . -  I  1 5  »  >  GOTO  40 
OK  =  .FALSE. 

TYPE  “LONGITUDE  INVALID,  MUST  BE  -90.  TO  -115, 

IF  ( . NOT. OK)  RETURN 

LASTME  =  TIME 

RETURN 

END 


,  TRY  AGAIN* 
IN" 

TRY  AGAIN'* 
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APPENDIX  B 


1.  Listing  of  job  control  stream,  JTCROTATE 

2.  Flow  chart  for  Program  TCROTAT 

3.  Listing  of  TCROTAT 

4.  Flow  chart  for  Subroutine  IMGPLOT 

5.  Listing  of  IMGPLOT 

6.  Flow  chart  for  Subroutine  ROTATE 

7.  Listing  of  ROTATE 


B-l 


3TCROTATE,  ID=UL 


VNULR ; T2000 ; HD 1 / EC215 > STSPC/SA  .  A6 
RFE (0) 

REMIND  (OUTPUT) 

ATTACH  /  ROTAT  ,  TCROT ATBIN ID=UL  . 
LIBRARY (KFNWCL IB) 

MAP (OFF) 

R E Q U E S T  <  I N  /  H D  ;  L  ,  N 0 R I N G  ,  V 3 N -71395) 

SKI  PF  ( I N  /ll#  1 7  /  B)*^J - 

COPYBF  (IN /TAPE-4) 

UNLOAD (IN) 

REWIND (ROTAT) 

RFE  (215) 

0N3W ( 1 ) 

ONSH (2) 

ROTAT. 

RFE (0) 

REWIND (TAPE 3) 

R  E  0  U  EST(OUT/HD/L/R  I N  6  V  3  N =71393)  <C 

SK I  PF  (OUT  /  8_/ 17  /  B )  - - - 

C  0  P  Y  B  F  ( T  A  P  E  3  /  0  U  T ) 

RETURN (OUT) 

*EOR 

222  269  --71 . 63  - 

it  EOR 
*EOF 


SUE  UHRICH 


VERIFY  VSN 

VERIFY  NUMBER  OF  FILES 
TO  SKIP 


VERIFY  VSN 

VERIFY  NUMBER  OF  FILES 
TO  SKIP 


INPUT  CARD 
In  this  case: 

ICNTR  =  222 
JCNTR  =  269 
ANGLE  =  -71.63° 
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FLOW  CHART 
FOR  PROGRAM  TCROTAT 


^  tcrotat  j 


V 


READ  SWITCHES 
READ  INPUT 
(SPIRAL  CENTER  & 

ANGLE) 


0 


1 


WRITE 

IMAGE 

ROTATED 

TO  A  NEW 

FILE 

\ 

r 

^  STOP  (9000^ 
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ooo  r~>oo  ooooooo 


PROGRAM  TCROTAT(INPUT«200B,OUTPUT«200B,TAPE3*0,TAPE4»0) 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THIS  PROGRAM  PROCESSES  SATELITF  IMAGES  THAT  ARE  ON  TAPE. 

THE  IMAGE  IS  READ  IN  AND  CENTERED  IN  THE  tyjINDOWt  ACCORDING  TO 
THE  INPUT  SPIRAL  CENTER.  THE  CENTERED  IMAGE  IS  THEN  ROTATED 
ABOUT  THIS  CENTER.  LASTLY,  THE  CENTERED  ROTATED  IMAGE  IS  WRITTEN 
BACK  TO  TAPE. 

INPUT: 

IMAGE  -  A  SATELITE  IMAGE  ON  FILE,  TAPEA  n 

ICNTR, JCNTR  -  PIXEL, RASTER  LOCATION  OF  SPIRAL  CENTER  OF  STORM. 

ANGLE  -  ANGLE  TO  ROTATE  IMAGE  COUNTER-CLOCKWISE  (DEGREES) 

°UTRIMAGE  -  AN  INTEGER  FILE  CONTAINING  THE  CENTERED  IMAGE 
PROCESSES : 

READ  IMAGE,  ANGLE,  AND  CENTER  LOCATION. 

CENTERS  IMAGE  TO  256,256  PIXEL  AND  RASTER  LOCATION 
ROTATES  IMAGE 

PRINT  ORIGINAL  AND/OR  IMAGE  ON  LINE-PRINTER 


SWITCH  SETTINGS: 
SWITCH  1  ON  = 
SWITCH  1  OFF  = 
SWITCH  2  ON  » 
SWITCH  2  OFF  * 


PRINT  OUT  ORIGINAL  IMAGE 
NO  PRINT  OUT 
PRINT  OUT  ROTATED  IMAGE 
DONtT  PRINT  ROTATED  IMAGE 


INPUT  FORMAT 
COLUMNS 
1-3 
5-7 
9-15 


PARAMETER 

X-VALUE  OF  CENTER,  LEFT-MOST  POINT  IS  1 
Y-VALUE  OF  CENTER,  TOP-MOST  POINT  IS  1 
ANGLE  TO  ROTATE  IMAGE  (DEGREES,  F7.2) 


CDMCONTAINSKTHEEFIRST  WORD  ADDRESSES  (FWA)  OF  EACH  RECORD  IN  THE 
IMAGE  FILES,  WHICH  ARE  STORED  IN  EXTENDED  CORE  STORAGE. 

I  EC  1  -  THE  INPUT  ORIGINAL  IMAGE 

I  EC  2  -  THE  ROTATED  IMAGE  DATA  (NO  HEADER  OR  D OC UME N A T I  ON ) 

INTEGER  BLOCK(IOO),  HE  ADER ( 69 ) 

INTEGER  IRASTER(512) 

COMMON  /ECBLK/  IECK515),  IEC2(512) 

- INITIALIZE  AND  READ  SWITCH  SETTINGS 

CALL  SSWTCHd,  ISW1  ) 

CALL  SSWTCH( 2, ISW2 ) 

- READ  INPUT  AND  CHECK  VALIDITY 

READ  111,  ICNTR, JCNTR, ANGLE 
IF  ( EOF ( 5L INPUT ) )  9000,100 

111  FORMAT( 13,  IX, I3,1X,F7.2)  ,  „„  „ 

100  IF  ( ICNTR. GT. 150  .AND.  I CN TR . L T . 362  )  GOTO  150 
PRINT  333, ICNTR, JCNTR, ANGL E 
333  FORMAT ( *  ERPOR  IN  INPUT:*/ 

1  5X,  ♦CENTER1'  *  ,  I  3  ,  *  ,  *  ,  I  3,  5X  ,  *  AN  GL  E  **  ,  F  7 . 2  ) 

GOTO  9000 

150  IF  ( JCNTR. GT .1 50  .AND.  JCNTR . LT .  362  )  GOTO  175 
PRINT  333, ICNTR, JCNTR, ANGL E 
GOTO  9000 

175  IF  ( ANGL E • L E . 1 80. 0  .AND.  A NGL E • G E •  —  1 80 . 0  )  GOTO  200 
PRINT  333, ICNTR, JCNTR, ANGLE 
200  PRINT  AAA,  ICNTR, J CNTR, ANG LE 
AAA  FORMAT ( *1  INPUT  DATA  CARD:*/ 

*2  *  CENTER  LOCATED  AT  *,I3,*,*,I3/ 

+  *  ANGLE  TO  ROTATE  IS*,F7.2,*  DEGREES*) 


C 

C 

C 

C 

c 

c 

c 


PEAD  IMAGE  FROM  TAPE,  CENTERING  IMAGE  AS  READ  IN.  THE  IMAGE  IS 
515  RECORD  BLOCKS  LONG,  AND  EACH  RECORD  IS  256  16-BIT  WORDS. 

THE  HEADER  RECORD  IS  READ  FIRST  AND  STORED  LAST.  THE  DATA  IS  ON 
RECORDS  2  TO  513  AND  IS  STORED  IN  RASTER  1  TO  512.  THE  DOCUMENATIC 
RECORDS  ARE  ON  514  AND  515  AND  ARE  STORED  AFTER  THE  DATA  RECOPDS. 
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c 

ID  =  2 56- I C NT R 
JD  =  256-JCNTR+l 
IBYT  -  I  A  B  S (  ID  )  +1 
I L  DC  =0 
C 

C  IF  JD  POSITIVE  MOVE  IMAGE  DOWN,  IF  NFGATIVE  MOVE  UP 
C  IF  ID  POSITIVE  MOVE  IMAGE  RIGHT,  IF  NEGATIVE  MOVE  LEFT 
C  FIRST  ZERO  OUT  OF  RANGE  RASTERS 
CALL  Z I L  CH ( BLOCK, 100) 

IF  (JD.GT.l)  GOTO  10 
IF  (JD.EQ.l)  GOTO  15 
1ST  *  512+JD 
DO  2  I”  1ST, 512 

I E  C 1  ( I ) -69* (  1-1 ) 

2  CALL  WR ITEC ( BLOCK, IEC1 ( I ) , 69 ) 

GOTO  15 

10  IED  =  JD-1  $  ILOC-JD-1 
DC  12  I  -  1 , I E  D 

IEC1 ( I ) =  6  9  *  ( 1-1 ) 

12  CALL  WRITEC ( BLOCK, IEC1 (I ),69) 

C  READ  IN  IMAGE 
15  REWIND  A 

BUFFER  IN  (9,1)  ( H E A DE R ( 1 ) , H E A DE R ( 69 ) ) 

IF  (UNIT ( 4  )  )  20,110,110 
20  IEC1 ( 515) =69*514 

CALL  WRITEC(HEADER,IEC1(515),69) 

DO  50  N-2, 515 

CALL  ZILCH( BLOCK  ,100) 

BUFFER  IN  (4,1)  (  B L OC K ( 1 ) , B L OC K ( 6 9  )  ) 

IF  ( UN  I T ( 4  )  )  25,110,110 
25  IF  ( N. GT. 513  )  GOTO  30 
IF  ( JD. LT. 1 )  GOTO  45 
I L  OC  -  ILDC  +  1 
IF  ( ILOC.GE.513)  GOTO  50 
IECl(ILOC)  =  6  9* ( I L  OC- 1 ) 

IF  ( ID. EO.O)  GOTO  40 

IF  (ID.LT.O)  CALL  BR PK ( B LOCK , 8 , I B YT , I R AS T E R, 8 , 1 , 5 1 2 ) 
IF  (ID.GT.O)  CALL  B R P K ( B L OC K , 8 , 1 ,  I R A S T E R , 8 ,  I  B YT , 5 1 2 ) 
CALL  WRITEC (  IRASTER, IEC1 ( ILOC ), 69) 

CALL  ZILCH(IRASTER,100) 

GOTO  45 

30  IECl(N-l)  -  69* ( N— 2 ) 

40  CALL  WRITEC(  BLOCK,  IECKN-1  ),  69) 

45  JD  *  JD+1 

50  CONTINUE 
GOTO  900 
C  ERROR  EXIT 
110  PRINT  777 
GOTO  10000 

FORMATS 

555  FORMAT ( *  ERROR  -  PARITY,  CONTINUING  EXECUTION*) 

777  FORMAT ( *  ERROR  -  EOF*) 


- PLOT  IMAGE  IF  SWITCHES  1  AND/OR  2  IS  ON,  ROTATE  IMAGE 

900  IF  (ISW1.E0.1)  CALL  IMGPLOT(l) 

CALL  ROTATE ( ANGLE  ), RETURNS (9000 ) 

IF  (ISW2.EQ.1)  CALL  IMGPL0T(2) 


WRITE  ROTATED  IMAGE  BACK  TO  TAPE.  (IN  FORMAT  READIBLE  BY  D.G.) 

- - WRITE  HEADER  BLOCK  TO  TAPE3 

N  =  1 

CALL  READEC(HEADER,IEC1(515),69) 

BUFFER  OUT  (3,1)  ( HE ADE R ( 1 ) , HE  A D E R ( 6 9  )  ) 

IF  ( UN  I T (  3  )  )  2020,  2110,2110 
C 

C - WRITE  DATA,  512  BLOCKS  LONG. 

2020  DO  2030  N=2,513 
LOC  -  N  —  1 

CALL  READEC (BLOCK, I E  C  2 ( LOC), 69) 
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BUFFER  OUT  (3,1)  ( B LOCK ( 1 ), BLOCK ( 69  )  ) 

IF  (UNIT ( 3 ) )  2030,2110, 2110 

2030  CONTINUE 

C-  _  -WRITE  DOCUMENTATION  BLOCKS  FROM  ORIGINAL  IMAGE 
DO  2040  N * 51 4 , 515 

L  OC  -  N-l  , 

CALL  READEC(BL0CK,IEC1(L0C),69) 

BUFFER  OUT  (3,1)  ( B L OC K ( 1 ) , B L OC K ( 69  )  ) 

IF  (UNIT ( 3 ) )  2040,2110,2110 
2040  CONTINUE 
ENDFILE  3 
GOTO  10000 

- ERROR  EXIT 


2110  PRINT  999, N 
999  F  0  RM  AT ( *  ERROR 
$  *  READ  1 

9000  PRINT  9999 
9999  FORMAT ( *  ERROR, 
10000  CONTINUE 
STOP 
END 


-  PARITY  OR  EOF  IN  WRITING  TAPE  -  NO. 
*,  I  5) 

STOPPING  EXECUTION*) 


BLOCKS*, 
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IMGPLOT 


Y 

PRINT  LEGEND 
PRINT  TITLE 


t _ 

READ  INTO  CM 
ALL  128  RASTERS 


FOR  EACH  RASTER 
UNPACK  DATA 
SELECT  PRINT 
CHARACTERS 
PRINT  RASTER 


RETURN 


non  o  ooo  oooo  ooooooooooooooo 


SUBROUTINE  I  MG  P L  OT  (ITYP) 


THIS  SUBROUTINE  PLOTS  THE  IMAGE  ON  THE  LINE-PRINTER 
INPUT : 

IMAGE  -  AN  INTEGER  FILE  CONTAINING  THE  INTENSITY  VALUE  OF 
EACH  PIXEL. 

ECBLK  -  COMMON  AREA  WITH  FW  A  OF  IMAGE  RASTERS 
IEC1  -  ORIGINAL  IMAGE 
IEC2  -  ROTATED  IMAGE 
ITYP  -  1  -  ORIGINAL  IMAGE  '  ' 

2  *  ROTATED  IMAGE 

OUTPUT: 

A  128X128  LINE-PRINTER  PLOT  OF  CENTRAL  PORTION  OF  IMAGE 

DIMENSION  IMAGE(128),  ICHARU6),  ITYPE<2) 

INTEGER  RASTFR(128),  BLOCK(88B2) 

COMMON  /  EC  B  L  K  /  IECK515),  IEC2(51 2) 

DATA  ICHAR/1H  , 1H . , 1 H J ,  1 H-, 1 H  +  , 1 H* , 1H1 , 1H 2 , 

S  1H3,1H4,1H5,1H6,1H7,1H8,1H9,1H=/ 

DATA  ITYPE/lOH  ORIGINAL  ,1CH  ROTATED  / 


- PRINT  TITLE  AND  LEGEND 

PRINT  111,  ITYPE(ITYP) 

111  FORMAT ( 1HO,  / 

1  59X,14H t*************  I 

2  59X, 1H/,12X,1H*/ 

3  59X,2H*  ,A10,2H  * / 

4  5  9  X  , 1  4  H  *  IMAGE  PLOT  */ 

5  59X, 1H*, 12X,1H*/ 

6  59X, 14H**** **********  ) 

PRINT  222 

222  FORMAT  (1HO/ 

1  53X,*LEGEND:*// 

2  53X»*INTENSITY  CHARACTER*/ 

3  5  3  X  ,  * -  - */) 

INTI  -  0 

INT2  *  15 
DO  700  I *1 » 16 

PRINT  333,  INTI,  INT2,  ICHAR(I) 

INTI  *  INTI  +  16 
INT2  *  I NT2  +  16 
700  CONTINUE 

333  FORMAT (53X, 14,*-*, I4,8X*A1) 

- ITYP  *  1,  ORIGINAL  IMAGE  -  ITYP  -  2,  ROTATED  IMAGE 

LW  «  -68 

PRINT  666,  ITYPE (  ITYP) 

666  FORMAT ( 1H1 , *FI LE  0 F* , A1 0 , * L T S T E D* ) 

IF  (ITYP.EQ.l)  I  EC  *1 EC1 ( 193  ) 

IF  (ITYP.E0.2)  IEC-IEC2( 193) 

CALL  READEC (BLOCK,  IEC,  8832  ) 

DO  1100  J  *  1 , 1 2  8 
LW  *  LW  +69 

- UNPACK  128  PIXELS  IN  A  RASTER  STARTING  WITH  THE  193RD  PIXEL. 

CALL  BRPK (BLOCK (LW),8,193, IMAGE, 60, 1, 128) 

K  =  0 

DO  1000  1*1,128 
K  *  K  +  1 
INT  *  IMAGE (  I) 

INT  -  INT  / 1 6  +1 

IF  (INT.LT.l  .OR.  INT.GT.16)  PRINT  888, INT 
IF  (INT.LT.l  .OR.  INT.GT.16)  INT  =  16 
RASTER(K)  ■  ICHAR(INT) 

1000  CONTINUE 

PRINT  777,  (RASTER(L),L*1,128) 

1100  CONTINUE 
777  FORMATdH  ,  1 2 8  (  A 1  )  ) 

888  FORM  AT ( *  ERROR,  INT**, 15) 

- ALL  DONE 

RETUPN 


B-8 


FLOW  CHART 

FOR  SUBROUTINE  ROTATE 


B-9 


non  or>o  ooooooooooooo 


SUBROUTINE  ROTATE ( ANGLE ) *PETURNS ( XXX) 

THIS  PROGRAM  ROTATES  THE  IMAGE  OF  SIZE  (M  X  N)  ABOUT  THE  CENTER 
BY  A  GIVEN  ANGLE  (WHERE  -1 8 0 . 1 < ANGL E < 1 80 . 1 ) 

INPUT : 

IMAGE  *  IMAGE  ON  MASS  STORAGE  FILE*  TAPE1 . 

ANGLE  =  THE  ANGLE  TO  ROTATE  THE  IMAGE  (COUNTERCLOCKWISE) 
ECBLK  =  COMMON  BLOCK  CONTAINING  FW A  OF  EACH  RASTER  IN  E.C.S. 

OUTPUT  : 

RIMAGE  *  ROTATED  IMAGE  ON  MASS  STORAGE  FILE*  TAPE2. 


INTEGER  IM  AGE ( 1 024  )*  RIMAGE(69)*  RASTER(130) 
INTEGER  R (512) 

COMMON  /ECBLK/  IECK515)*  I  EC2  t  512  ) 

- ROTATE  AXIS  ABOUT  C EN TE R (  2 5 5 . >  2 5 5  .  ) 

IF  ( ANGLE. EO. 0.0)  RETURN  XXX 
P  DR  D  =  255.0  *  255.0 
ARAD  =  ANGLE  *  0.01745329 
SINANG  =  SIN(ARAD) 

COSANG  -  COS(APAD) 

DO  200  J  *  1 *  512 

CALL  ZILCH(R*512) 

I  EC  2 ( J ) =69* ( 514  +  J  ) 

Y  =  J-l 

YOLD  =  255. -Y 

X  2  X  2  =  R  DP  D  -  (  YO  L  D  )  *  (  Y  OL  D  ) 

IF  (X2X2.LT.0.)  GOTO  100 
X 2  =  SORT  (X2X2)  +  255. 

12  =  X  2  +  1.0 

II  «  511.5  -  X 2 

IF  (  12. GT.  ID  GOTO  125 
100  CALL  WR ITEC ( R* I EC2 ( J)*69) 

GOTO  200 

125  DO  150  1=11*12 

-  XOLD  =  1-256 

- CALCULATE  NEW  POINT  XPRIME*  YPRIME 

XPRIME  =  ( XOLD+COS ANG)  +  ( YOL D* S I N ANG ) 
YPRIME  =  ( YO LD*COS ANG )  -  ( XO L D* S I N AN G ) 

IX  =  X  =  XPRIME  +  255. 

IY  =  Y  *  255.  -  YPRIME 

IF  (IX.LT.l  .OR.  IX.GE.512)  GOTO  150 

IF  (IY.LT.l  .OR.  IY.GE.512)  GOTO  150 


— INTERPOLATE 

INTENSITIES 

CALL 

READEC ( RASTER*  I 

CALL 

BRPK ( RA  STER* 8, 1 

CALL 

BRPK ( R  A  STER (70) 

PO  = 

I M  A  G  E ( IX) 

PI  » 

IMAGE ( IX+1) 

P2  = 

IMAGE ( I X+512 ) 

P3  = 

IMAGE ( I X+ 51 3) 

DX  * 

X  -  FLOAT(IX) 

DY  = 

Y  -  FLOAT(IY) 

RINT 

»  ( P 1-P  0 ) *  D  X  + 

S 

( P3+P0-P2-P 

R  (  I ) 

=  RINT+.5 

150 


( P  2  — PO ) *  DY  + 

1 ) *  DX*  DY  +  PO 


CONTINUE 

CALL  BRPK ( R* 60* 1* RI M AGE  *  8* 1* 51 2 ) 
CALL  WRITEC (RIMAGE* IEC2 (J ), 69) 
200  CONTINUE 
RETURN 
END 
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1.  Listing  of  job  control  stream,  ^INTERPOLATE 

2.  Flow  chart  for  Program  ITRPIMG 

3.  Listing  of  ITRPIMG 

4.  Flow  chart  for  Subroutine  READIMG 

5.  Listing  of  READIMG 

6.  Flow  chart  for  Subroutine  SUMSN 

7.  Listing  of  SUMSN 

8.  Flow  chart  for  Subroutine  ROTAVG 

9.  Listing  of  ROTAVG 

10.  Flow  chart  for  Subroutine  INTRP 

11.  Listing  of  INTRP 

12.  Flow  chart  for  Subroutine  PRT 

13.  Listing  of  PRT 

14.  Flow  chart  for  Subroutine  RITIMG 

15.  Listing  of  RITIMG 
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JINTERPOLATE,  ID  =  UL 


VNULI  *  T 4 0 0 0 *  EC320*  HD1  *  *  STSF'C*  S6.  46 
RFE? 0. 

REWIND* OUTPUT 
REWIND* OUTPUT. 

ATTACH? I  TRP* INTERPOLATEBIN* ID-UL. 

L I BRARY  * ♦FNWCL I B . 

map*dff. 

REQUEST *  I N  ?  VSN=7 1 396 *  RD  *  L  *  NOR I N6 . 

SKIPF* IN* 18* 17*  B.  - - 

COPYBF* IN*  TRPE3. 

COPYBF *  IN*  TRPE4. 

REW I ND  *  TRPE3 *  TRPE4 . 

UNLOAD* IN. 

□NWS*  1 .  - - - - - 

RFE*  380. 

I TRP . 

RFE*  0. 

RETURN* I TRP. 

REWIND* TAPES. 

REQUEST  *  OUT  *  VSN=71397  *  HD  *  L  *  R I NG  .-^■ 

SK I PF  *  OUT  *jB. *  1 7  >  B .  - - - 

COPYBF* TAPES* OUT. 

RETURN* OUT. 

EXITjS. 

♦EOR 

150  - - - - - - 

♦EOR 

♦EOF 


VERIFY  VSN 

VERIFY  NUMBER  OF  FILES  TO  SKIP 


PRINT 


VERIFY  VSN 

VERIFY  NUMBER  OF  FILES  10  SKIP 


INPUT  CARD*  RADIUS  OF  CROSS- 

CORRELATION  CIRCLE 
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ITRPIMG 


_ I _ 

CALCULATE 

ROTATION 

ANGLE 
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PROGRAM  I TRP I  MG  <  I NPUT =£  0  QB  *  OUTPUT =2  0  0 B  *  T  APE3= 0  *  TAPE4=0*  TAPE5=0> 

C 

C  THIS  PROGRAM  CRERTES  AN  I WAGE?  WHICH  IS  INTERPOLATED  FROM  TWO 
C  SAT ILL I TE  IMAGES  READ  FROM  TAPE.  FROM  A  RADIUS  OF  AM  AMMULUS 
C  READ  IN*  CROSS  CORRELATION  IS  PERFORMED  TO  DETERMINE  THE  OPTIMUM 
C  ANGLE  TO  ROTATE  THE  IMAGES  INPUT  TO  CREATE  THE  NEW  IMAGE. 

C 

C  INPUT  WILL  BE  TWO  IMAGES*  AND  THE  RADIUS  SPECIFIED  ON  THE  INPUT  CARD. 

C 

C  OUTPUT:  WILL  BE  THE  INTERPOLATED  IMAGE  AND  ANY  CALCULATIONS  MADE 
C  DURING  THE  RUN  OF  THIS  PROGRAM. 

C 

C  SWITCH  1  ON:  PRINT  OUT  BEFORE  AND  AFTER  IMAGE*  AND  ALSO 
C  PRINT  CENTRAL  PORTION  OF  INTERPOLATED  IMAGE  (.128X1  £8)  - 

n 

COMMON  /ECBLK/  I  EC  1  <5 1 5>  *  I  ECS  <5 1 5>  *  I  ECS  <5 1 2> 

COMMON  /  INTENSE  IMG1  (1602)  *  I  MG 2  <  1 6 0£> 


C 


c 

c 


c 

c 

c 


CALL  SSWTCHC1*  IS1':> 

PRINT  1 

1  FORMAT  <1H0*  T£5*  29H+-  PROGRAM  INTERPOLATE  ♦// 

%  IX*  T1 1  *  50HTHIS  PROGRAM  WILL  INTERPOLATE  AN  IMAGE  BETWEEN  TWO 

$  *3H  OTHERS. /> 

— READ  INPUT  AND  VERIFY 
READ  2* I RAD I US 

2  FORMAT <1 3> 

PRINT  3*  I  RADIUS 

3  FORMAT <1H0* ♦RADIUS  EQUALS  ♦» I 3> 

IF  < I RAD I US. GE. 1  .AND.  I RADIUS. LE. £55)  GOTO  25 
PRINT  4 

4  FORMAT C1H0* 50HRADIUS  MUST  BE  GREATER  THAN  ZERO  AND  LESS  THAN  £55/ 

+  1H  *  21HPR0GRAM  STOPPING  NOW. > 

GOTO  90 


-  -READ  IN  IMAGES  AND  PRINT  IF  SWITCH  1  IS  ON 

25  CALL  READ I MG* RETURNS <90> 

IF  C1S1.EQ.2)  GOTO  30 
CALL  PRT  <1> 

CALL  PRT  <2> 


C  CONVERT  X *  Y  LOCATIONS  ON  CIRCLE  TO  R* THETA  COORDINATES. 

C  CORRESPONDING  INTENSITIES  ARE  STORED  IN  THE  COMMON  BLOCK*  INTENS. 

C 

30  RAD  =  I RADIUS 
DTHETA  =  1 . /RAD 
NOPTS  =  2.  3.14159265  ♦  PAD 

PRINT  5*  NOPTS/ 

5  FORMAT  OONUMBER  OF  POINTS  ON  CIRCLE  TO  CONVERT  =♦* 15) 

DO  40  1=1* NOPTS 

THETA  =  C I -1> ♦DTHETA 
X  =  RAD  ♦  COS  <  THETA)  +  255.0 
Y  =  255.0  -  CRAD  ♦  SIN  C THETA  -  ) 
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c-  -  -STQRE  INTENSITY  IN  ARRAY  BY  I  EQIJRL  THETR  STEP  SIZE 

C - 1  URGE  1 

CRLL  I NTRP  OjXjYjI  NT) 

IMG1 <I)=INT 

C - IMAGE  2 

CRLL  I  MTRP  C2 1.  X  j  Y I  MT) 

IMG2  Cl) =IMT 
40  CONTINUE 
C 

C  CALCULATE  CROSS-CORRELATION.  DONE  IN  TWO  STEPS. 

C  NINC  =  INCREMENT  OF  STEP  SIZE 

C  NMAX  =  VALUE  OF  N-  RT  MAX  CORRELATION 

C 

NMAX  =  0 
NINC  =  5 

CRLL  SUMSN CNOPTS • N I NC » NMAX) 

NINC  =  1 

CRLL  SUMSN  CNOPTS > N I NC  t NMAX) 

C 

C - COMPUTE  RNGLE  AND  ROTATE  IMAGES 

ANGLE  =  NMAX  ♦  DTHETR'S. 

A  =  ANGLE--'.  01745329 
PRINT  3j  ANGLE  j  A 

S  FORMAT OONOW  ROTATING  IMAGE  A  BY  ♦  jF7.4>*  RADIANS  OR  ♦jF7.£j 
$  *  DEGREES*) 

IF  '"ANGLE.  EQ.  0.  >  GOTO  90 
CALL  ROTAVG < ANGLE) 

IF  Cl S 1 .  EQ .  1 )  CALL  PRT  C3) 

C  WRITE  NEW  IMAGE  TD  TAPE 
CALL  RITIMG 
90  CONTINUE 
100  STOP 
END 
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n  o  o  n  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  fit 


SUBROUTINE  RE RBI MG .RETURNS (XXX) 

THIS  SUBROUTINE  READS  IN  TWO  I URGE  FILES.  THE  BEFORE  INRGE  IS  ON 
TRPE3  RNB  THE  AFTER  INRGE  IS  ON  TRPE4. 


FDRNRT  OF  EACH  INRGE  FILE  <RECDRD 
RECORD  NO.  CONTENTS 

1  HERDER 

2-513  DRTR 

514-515  BDCURENRTIDN 


*  256  16-BIT  WORDS> 


FORNRT  OF  OUTPUT 
RECORD  ND. 

1-512 

513-514 

515 


FILES  (REC0RB*69  WORDS.  >: 

CONTENTS 

DRTR  <512  8-BIT  PIXELS  PER  RECORD) 
DDCUNENTRT I ON  (UNCHANGED) 

HERDER  (UNCHANGED  RECORD  FORNRT) 


D I NENS I  ON  I NRGE (69) . I HERDER  <1 6  0) . I WK  <1 6 0> 
CONNON  ✓ ECBLKs  IEC<1030)>  IEC3C512) 

- READ  TWO  INRGE  FILES 


REWIND  4 
REWIND  3 

NW=1  *  I FILE=3 
5  N=1 

STORE  HERDER  RECORD  LAST 

BUFFER  IN  (IFILE. 1)  (I NRGE  <1 > . I NRGE (69> > 

IF  (UNIT (IFILE) >  10.110.110 
10  I  EC  (NW+514)  =  69^<NW+513> 

CALL  WR I TEC  < I NRGE » I EC (NW+5 1 4) . 69) 

CALL  BRPK (INRGE. 8. 1. I HERDER. 12.1.80) 

CALL  RSC I DPC < I HERDER. 1 6 . I WK . 1 6 . LDPC . I ERR) 
IF  (I ERR. EQ. 1)  PRINT  1 1 1 . LDPC 
PRINT  222. (IWK(L) .L=l. LDPC) 

DO  100  N=2. 515 

I EC (NW) =69* (NW- 1 ) 

BUFFER  IN  (IF I LE . 1 )  (I NRGE (I) . I NRGE  <69) ) 
IF  (UNIT  (IFILE))  50.110.110 

- STORE  DRTR  AND  DDCUNENTRT ION  RECORDS 

50  CALL  WRITEC (INRGE. IEC (NW) . 69) 

NW=NW+T 
100  CONTINUE 

IFILE=IFILE*1 

NW=NW+1 

IF  (IFILE. EQ. 4)  GOTO  5 
RETURN 


- ERROR  EXIT 

110  PRINT  333.  IFILE.N 
RETURN  XXX 

- FORNATS 


111  FORNRT  O  ERROR  IN  ASCII  TO  DISPLAY  CODE  CONVERSION.  LDPC  =*.16) 

222  FORNRT  <♦*  INRGE  FILE  HERDER  IS  *.8R10) 

333  FORNRT  <*  ERRDR  —  PARITY  OR  EOF  IN  READING  —  STOPPING  EXECUTION*’-" 
+  ♦  FILE  NUNBER*’.  12.*-  NUNBER  OF  BLOCKS  READ*-.  15) 

END 
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FLOW  CHART 

FOR  SUBROUTINE  SUMSN 
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o  o  o  o  o  o 


SUBROUTINE  SUNSN CNDPTS >  N I NC »  NMRX.> 

THIS  SUBROUTINE  CRLCULRTES  THE  CRQSS-CDRRELRT I ON?  S CN> r  WHICH  IS 
(THE  SUN  OF>  IMG!  <I>  _TIMES  ING2a*N>  FOR  N  E0URL  TO  +✓-  NINOS. 

INPUT:  NINC*  INCRENENT  DF  PIXELS 

NMRX*  ZERO  DR  THE  NINC  WHICH  GRVE  THE  LRRGEST  CORRELRTION 
NOPTS,  NUMBER  OF  POINTS  IN  ERCH  CONVERTED  I URGE  RRRRY 
INTENS*  COMMON  RRER  CDNTRINING  THE  INTENSITIES  OF  THE  CONVERTED 
C  PIXELS  ON  THE  RNNULUS. 

C  OUTPUT:  NMRX  THE  HIGHEST  SN  VRLUE  RT  THIS  N 
C 

DINENSION  SNC1 1 > 

CONNON  -'INTENSE  INSta602>*  IMG2C1602> 

C 

10  SNNRX  =  0. 

NN  =  N  =  — 5*NINC  +  NNRX 
PRINT  HI?  NINC 

111  FDRNRT <  1  HO  * 35HN  VS.  S (N>  *  WHICH  INCREMENT  OF  N  IS? 13) 

PRINT  112 

112  FORMRT C1H0* ♦  N  ♦*  T20*  +$■  <N>  ♦»"> 

C 

DO  100  IN=1*11 
SN  <IN>  =  0. 

DD  50  J=1 j  NDPTS 
K  =  J+N 

IF  CK.LT.l*  K=K>NOPTS 

IF  <K. 6T . NOPTS>  K=K-NDPTS 

SN<IN>  =  IMG1  <J>  ♦'ING2  <K>  +  SN<IN> 

50  CONTINUE 

PRINT  114;.  N*  SN<IN> 

IF  CSN < I N> . LT. SNNRX>  GOTO  30 
SNNRX  =  SNCIN> 

NNRX  =  N 

30  N  =  N+NINC 

100  CONTINUE 

114  FORMRT  < I H  ? 1 3  ? F22 . 2> 

PRINT  1 15* NNRX 

115  FORMRT  <1  HO* +-  N  NRXINUN  IS  *■>  I3> 

IF  CNN.  EQ.  NMRX  .OR.  IRBS  <NN> . EQ. NNRX>  GOTO  10 

RETURN 

END 


C-1  0 


FLOW  CHART 

FOR  SUBROUTINE  ROTAVG 
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o  n  o  o  o  o  o  o  o  o  o 


SUBROUT IME  ROTAVG  < AMGLE> 

,0  i 

m  THIS  *  SUBROUTINE  ROTATES  IMAGE!  OMCIWISE  AMD  I  URGES  COUHTER- 
%;t  CLOCKWISE  THE  SAME  AMOUNT.  THE  IMAGES  ARE  THEN  AVERAGE B  TO  GET 

'*'y 'time.  ■TllP*  yrtl  TiuVf-nnm  A-Virm  irwu-r> 


5Ff  THE  HEW  INTERPOLATED  IMAGE 

'i'i  *  >1'^"  ”  ■  f  \ 


INPUT:  AHGLE  -  AMOUNT  IH  RADIAHS  TO  ROTATE  THE  IMAGES 

ECBLK  -  FIRST  WORD  ADDRESSES  OF  EACH  RASTER  OF  THE  IMAGE 

OUTPUT:  THE  HEW  IHTERPOLATEB  IMAGE 


DIMENSION  IMG3<69>.  NEW<512>  .  INT<2>.  COSAHG  <2>  ,  SIHAMG  <2) 
COMMON  /ECBLK^  IEC1<515>.  IEC2<515>.  IEC3<512> 

- IMITIRLIZE  VARIABLES 


S I HAHG  <  1  >  =  SIN  < ANGLE) 
COSAHG a >  =  COS  < ANGLE) 

S'  I  HAHG  <2)  =  SIH  <— ANGLE) 
COSAHG <2>  =  COS  <-ANGLE) 
RR  =  255.  255. 

- ROTATE  IMAGES 


DO  200  J  =  1 . 512 

C-  -  -CALCULATE  RASTER  DF  MEW  IMAGE;.  JUST  DO  A  255.0  RADIUS  CIRCLE 
IEC3  <J>  =69*<J+I  029) 

CALL  ZILCH <NEW. 512) 

Y  *  J  -  1 

IF  <Y. EQ. 511.)  GOTO  50 

YOLD  =  255.  0  -  Y 

X2X2  =  RR  -  <Y0LD>  ♦■CYOLD) 

IF  <X2X2. LT. 0. )  GOTO  50 
X2  =  SORT  <X2X2)  +  255.0 
12  =  X2+1. 

II  =  511.5  -  X2 
IF  <12. GT. I 1>  GOTO  100 
50  CALL  WRITEC  <NEW. IEC3  <  J)  j  69> 

GOTO  200 

C - IHTERPOLATE  DHE  RASTER  DF  EACH  IMAGE 

100  DO  150  1=11. 12 

XOLD  =  I  -  256 
DO  125  IM=T . 2 

XPRIME  =  <XOLD*COSANG<IN>)  +  < YOLD*S I  HAHG  <  I M)  > 

YPRIME  =  <YOLD*COSANG  <IN> )  -  <XOLD*SINANG <IM> ) 

X  =  XPRIME  +  255.0 
Y  =  255.0  -  YPRIME 

IF  (X.LT.1.0  .OR.  X.  GE.  512.  0)  GOTO  150 
IF  CY.LT.1.0  .OR.  X. GE. 512. 0>  GOTO  150 
CALL  I MTRP  < IN. X  ?  Y . I NT  <  I M>  > 

125  COHTIHUE 

C - AVERAGE  PIXELS 

IF  <INT<l6*EQ.  0)  IHT  <2> -0 
IF  <IHT  <2) . EQ. 0>  IHT  <  I )  =0 
HEW  <I>  ='  FLOAT  <INT  <I) +TNT  <2>>  ♦>  0.5  +.5 
150  COHTIHUE 

C - PACK  RASTER  <512  PIXELS  IHTO  63  WORDS)  AMD  STORE  RASTER. 

CALL  BRPK <HEW. 60. 1 . IMG3. 8. 1 .512) 

CALL  WR I TEC  <1 MG3 . 1 EC 3 < J> . 69) 

200  COHTIHUE 
RETURM 
EMD 
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o  o  o  o  o  o  ^  '  ji  o  o  o  o  o  o  o  o  o  o  o  o  n  o  o  o 


SUBROUTINE  I NTRP CIFILErX* Yr  INT> 

THIS  SUBROUTINE  DOES  A  BILINE HR  INTERPOLATION  TO  GET  THE  INTENSITY 
VALUE  AT  POINT <XyY> 

INPUT? . 

I FILE  -  FILE  TO  READ  RASTERS  FRON 
Xy  Y  -  LOCATION  OF  POINT 

FRON  TOP-LEFT  CORNER  OF  INAGEy  BEING  (M.) 

ECBLK  -  CONNON  BLOCK  WITH  STORAGE  LOCATIONS  OF  EACH  RASTER 

OUTPUT • 

INT  -  INTERPOLATED  INTENSITY  VALUE  AT  <.X>Y) 

DINENSION  IRf<512>y  IR2<512)y  IDATAC138) 

CONNON  /ECBLK^  IECI <515)  y IEC2 <31 5> > IEC3 <512) 

- FIND  LOCATION  OF  INTENSITY  VALUES 

IY  =  I  FI  X  <Y) 

IX  =  I F I X  <X> 

IF  <IFILE.EQ.2>  GOTO  5 
CALL  RE ADEC < I DAT Ay  I EC 1 < I Y>  » 1 38) 

GOTO  10 

CALL  READEC  < I DATA* IEC2<IY>  y 138) 

0  CALL  BRPK <IDATAy8y 1 y IRlyGOy 1 y512) 

CALL  BRPK<IDATA<70>  y8y  ly  IR2y60y  ly512) 


- GET  INTENSITIES 

PO  =  IR1 <IX) 

PI  *  IR1 <IX+1> 

P2  =  IR2  <IX) 

P3  =  IR2<IX+1> 

- CALCULATE  INTENSITY  AT  POINT <XyY)  FRON  SURROUNDING  POINTS 

DX  =  X  -  FLOAT <IX) 

DY  *  Y  -  FLOAT  <IFIX  <Y)  ) 

XI  NT  =  (Pl-P  0>  ♦BX  +  <P2-P0)+DY  +•  <P3+P0-P2-P1>  *DX*DY  +  PO 

INT  =  XINT  +•  0.5 

RETURN 

END 
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o  o  o  o  o  o  o  o  o  o  o 


SUBROUT I ME  PRT<IHGFILE> 


THIS  SUBROUT I ME  PR I MTS  OUT  THE  CENTRAL  128X1 £8  PORTION  OF  RM  I  URGE. 

INPUT?.  IMGFILE  -  UNIT  NUMBER  OF  THE  FILE 

UNIT  1  -  BEFORE  II1BGE 
£  -  RFTER  I URGE 
3  -  INTERPOLATED  IMAGE 
ECBLK  -  STORAGE  LOCATION  OF  EACH  RASTER 

DIMENSION  IMAGE <128>  ?  IM<:i£8>?  ICHAR<16>?  NAME  <3>  ?  IFILE«.883£> 
COMMON  /ECBLK/  IEC1 <515> ? 1EC2 <515> ? IEC3 <512> 

DATA  I  CHAR/ 1 H  ?  1H.  ?  1H?  ?  1H-?  1H+j  1H*-?  1H1 ?  1H2? 

+  1H3?  1H4?  1  H5?  1H6?  147?  1H8?  1H9?  1Hs,' 

DATA  NRME/1  OHTHE  BEFORE;.  1  OH  THE  AFTER?  1  OH  THE  NEW/ 


PRINT  1?  NAME < I MGF I LE> 

1  FORMAT  <  1  HI  ?A10?+-  IMAGE  PRINTED** 

IF  <IMGFILE. EQ. 1>  IEC=IEC1 <193> 

IF  < IMGFILE. EQ.£>  IEC=IEC2 <193> 

IF  < I MGF ILE . EQ . 3>  IEC=IEC3<I93> 

CALL  READEC  <IFILE? IEC?  8832> 

L=T 

DD  400  J= 193? 320 

CALL  BRPK  < I FILE <L>  ?  3? 1 93 ? IMAGE?60?1? 1£8> 
DO  300  1=1  ?1 £8 
K  =  I 

I NT  =  IMAGE < I > 

I NT  -  I NT/ I 6  +  1 
IM<K>  =  I CHAR  < I NT> 

300  CONTINUE 

PRI  NT  £ ?  <  I M <K>  ?  K=1  ?  1 28> 

L=L+69 
400  CONTINUE 

2  FORMAT  <1 H  ? 1 28  <A1 > > 

RETURN 

END 
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017 


i_  ^  stjBRGUTIRE  tmtie*  ■>"' 

THIS  SUBR£UfTHE  'WRITES  THE  HEWLY  CREATED  IHfiGE  TO  THE  FILE?  TRPE5. 
€  IMPUTr- 

C  ECBLK  -  COHHOH  BLOCK  WITH  STORAGE  LOCRTIOHS  OF  ERCR  RRSTER 
C  IEC1  -  BEFORE  IHRGE 

C  IEC2  -  RFTER  IHRGE 

C  IECG  -  HEW  IHTERPOLHTED  IHRGE 

C  OUTPUT f  THE  HEW  IHTERFDLRTEB  IHRGE  WRITTEH  TO  TRPE  OR  DISK  FILE. 

C  THE  FQRHRT  DF  THE  FILE  WILL  BE: 

C  RECORD  HO.  CDHTERTS 


C 

C 

c 

c 

c 

c 

e 

c 


1  HERDER 

2-513  DRTR  (PIXEL  IHTEHSITIES> 

514  -  515  DOCUMEHTRTIOH 

THE  HERDER  RHD  DDCUHEHTRT I  OH  RECORDS  WILL  BE  COPIED  FROF1 
THE  BEFORE  IHRGE. 


D I  HERS  I  OH  I HRGE  (69>  -  I  HERDER  (69>  ?  I WK  (8> 

CDHHOH  /ECBLK/  IEC1 (515> *  IEC2(515>?  IEC3(512> 


DRTR  IWK/t  0HIHRGE4IHTE?  1  OHRPOLRTED4R?  1  0HFTER444444?  5+0/ 


C 

C - IHITIRLIZE 

C 

RSKt  =  1001 001 001 001 001 001  OB 
HSK2  -  01 001 001 001 001 001  001 B 
HSK3  =  001 001 001 001001 001 00B 


C 

C - HERDER  RECORD 

C 

CRLL  BRPK  ( I WK  ?  6  ?  1  ?  I WK  ?  9?  1  ?  3 0> 


I WK  ( 1 >  -  I WK  ( 1 > 

.OR. 

HSK 1 

I WK (2>  =  IWK (2> 

.  OR. 

HSK2 

IWK (3>  =  IWK (3> 

.OR. 

HSK3 

I  WK  (4>  =  I  WK  (4> 

.OR. 

HSK  I 

CRLL  BRPK  (IWK?9? 1 ? I HERDER?  8?  1  ?  30> 

CRLL  RERDEC  (IHRGE? IECt (515>  ?69> 

I HERDER (4>  -  SHIFT  (  IHRGE (1> . RHD. HRSK (56> ? +56> 
ICHRR  =  IHRGE  <I>  .  RHD.  .  HOT.  HRSK  (5G> 

ISCRTCH  =  I HRSE (2> . RHD. HRSK <56> .OR. ICHRR 
I HERDER (5>  =  SHIFT  (ISCRTCH:.  +56> 

CRLL  BRPK  < I HERDER?  8 *  1 ? I WK ? 1 2 1 1 ? 8  0> 

CRLL  RSCIDPCdWK?  16?  IHRGE?  16?LDPC?  IERR> 

IF  (I ERR. EQ. 1>  PRIHT  222?LDPC 
PRIHT  333?  (IHRGE <L>  ?L=1?LDPC> 

I HERDER <5>  *  IHRGE (2> 

BUFFER  OUT  (5 ? 1 >  (I HERDER ( 1 >  ? 1 HERDER (69> > 

IF  (UHIT  (5>  >  20?  110?  110 
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o  o  o  o  ct 


-  -v-^ITE  BATA  HMD  DDCUMEMTATIOM  RECORDS 

20  DO  30  N=2?515 
LOC  -  H-t 

IF  <L0C.GT.512>  GOTO  25 

CALL  READEC  <  I MAGE?  IEC3  <LDO  9  69> 

GOTO  27 

25  CALL  READEC < I MAGE? I EC 1  <LDC> ?  63> 

27  BUFFER  DUT  <5?  1  >  <  I  MAGE  <t>  ?  IMAGE  <69» 

IF  <UNIT<3>>  30.1 10i.H0 
,  30  COMTIMUE 

EMDFILE  5 
REWIMD  5 
RETURM 

- ERROR  EXIT 

C 

no  primt  Hurt 

111  FORMAT <♦=  ERROR  —  EOF  OR  PARITY  IM  RITIMG?  MO.  BLOCKS  =  ♦.IS 
222  FORMAT <♦*  ERROR  TDD  SMALL  IMAGE  ARRAY?  LDPC  =**?I5> 

333  FORMATS  MEW  IMAGE  HEADER  IS  *?  8A1 0> 


C 


RETURM 

EMD 
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Listing  of  Abbreviations  and  Definitions  used  in  this  report 
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ABBREVIATIONS  AND  DEFINITIONS 


CDC 

Control  Data  Corporation 

CM 

Central  Memory 

DG 

Data  General 

ECS 

Extended  Core  Storage 

FNOC 

Fleet  Numerical  Oceanography  Center 

GADHS 

GOES  Acquisition  and  Data  Handling  System 

GOES 

Geostationary  Operational  Environmental  Satellite 

HAL 

A  computer  at  FNOC  (a  CDC  6500) 

MCIDAS 

Man  Computer  Interactive  Data  Access  System 

NEPRF 

Naval  Environmental  Prediction  Research  Facility 

NTSC 

National  Television  System  Committee 

ODSI 

Ocean  Data  Systems,  Inc. 

PGP 

Programmable  Graphic  Processor,  part  of  the  Genisco 
GCT-3000  display  unit 

Pixel 

Picture  element,  the  smallest  addressable  element  of  an 
image 

Raster 

A  row  of  pixels 

RBG 

Red-Blue-Green,  the  basic  colors  from  which  digital 
images  are  built 

SPADS 

Satellite-data  Processing  and  Display  System 

SPC 

A  computer  at  FNOC  (a  CDC  CY1 75) 

VSN 

Volume  Serial  Number,  for  magnetic  tapes 
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